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Abstract 



'~t In the context of some bidimensionnal Navier-Stokes model, we exhibit a family of exact 

'""5 oscillating solutions {"Ueje defined on some strip [0, T] x R^ which does not depend on e €]0, 1]. 

iy-\ The exact solutions is described thanks to a complete expansions which reveal a boundary layer 

in time t — 0. The interactions of the various scales (1, 1/e and 1/e^) produce a macroscopic 

effect given by the addition of a diffusion. To justify the existence of {ue}e, we need to perform 

various Sobolev estimates that rely on a refined balance between the informations coming from 



■^^ the hyperbolic and parabolic parts of the equations. 

^ Mathematics subject classification (2000): 35-XX, 76N17, 76M45 
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B 

1 Introduction 

^ In Section 1, we introduce the underlying equations and the functional framework. Then, we state 

C*^ our main result. 

00 

O 1-1 The equations 

^^ The time and space variables are t e M_|_ and x := (a;i, 0:2) G M^. The state variables are the density 

^•^ p e M+ and the two components ui and U2 of the velocity of the fluid u :— *(ui, U2) G IR^- Given a 

,__! function u : M^ — > M^, note as usual: 

> d d 

• ^ div u ■- diui + d2U2 , di -.^ - — , 82 := t^ — . 

k>( ox I dx2 

^ In what follows, e G 10, 11 is a parameter approaching zero. Introduce the dissipation: 

Veu = * {Vlu , V'^.u) := p. e^ A^jU + A e^ V div u 
where /u, A G M^ are fixed. Let h be some smooth periodic function with mean zero: 

h:T — >M, T:=M/Z, /igC°°(T;M), f h{e) dO = . 

Jf 
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Consider the following oscillation which is polarized on the second component: 

F,{x) = 'iO,F^)ix) := £-2 '(o,^i^^eeh {e~^xi)) , e e]0,l] . 

Our starting point is the study of a model based on two-dimensional compressible isentropic equa- 
tions of the Navier-Stokes type, as can be found in [3, 12], forced by the source term F^ : 

\dtp + div{pu)^0, fuiui uiU2 

\dtipu)+divipu(g)u)+Vp''' ^ p{VeU~ F^) , \U1U2 U2U2 



where 7 the adiabatic constant is supposed to be larger than one. To obtain a quasi linear system 

^p^ with C := ^ 



having a symmetric form, it is classical [13] to introduce the state variable p := ^p'^ with C :— - — 



Then, we have to deal with: 

{dtp + u -^/p + C pdbfu — Q ^ , . 

dtu + u-Vu + CpVp^^VeU- F^. 

Observe that: 

V/{0,h,)-F,^0, h,{x) -.^ h(^^y V£e]0,l]. 

It follows that, for all e g]0,1], the oscillation *(0,0, /i^) satisfies Equation (1). 

Our aim is to consider the problem of the stability of such families of solutions. To this end, at 
the initial time i = 0, we modify *(0, 0, h^) by adding some perturbation. More precisely, we start 
with: 

{p,u\u^){0,x) ^ (0,0,/.)(^) + (e'^go.e,e*^<e,£"'<e)(5,y) (2) 

where {i^,M) E N^ with i^ large enough and M > 7/2 (retain that i/ :^ M), whereas: 

{qo,e,vl,,vl,){9,y) eH°°{TxR;R'), y.^^eR. 

One effect of the above perturbation is to introduce a dependence on X2 G M (or y £R). Despite 
the smallness of e" (and maybe e^^), when solving (l)-(2), we have to understand the interactions 
that occur between the very fast oscillations in the direction xi (with wavelength e"^) and the fast 
variations in the transversal direction X2 (with wavelength e). On this way, we are faced with 
questions about turbulence, in the spirit of models proposed in [4, 5, 7]. 

Another insight on the subject can be obtained by looking at (1) in the variables {6,y) G 
T X M. Then, we are faced with a hyperbolic-parabolic system implying some singular (in e G ]0, 1]) 
symmetric quasilinear part: 

dtp + e^"^ (u^ dep + e V? dyp) + C e"'^ p (dgu^ + e dyu"^) = , 

dtu^+e-^ lu^deu^+eu^dyU^)+Ce-^pdep^V}u, (3) 

and some viscosity which is degenerate on the density and becomes large when e — >■ 0: 



V.u 




p {degu^ + e^dyyu^) + Xe {degu^ + edeyU^) 
p {dggu^ + e^dyyV?) + Xe {edeyU^ + e'^dyyU^] 



In this article, we show that (for v large enough and M > 7/2) the oscillating Cauchy problem 
(l)-(2) is locally well posed in time. We prove (Theorem 1.4) the existence of a time T G K+ 
independent of e S]0, 1] with solutions {p^,ul,u^) = {e'^ qg,e^ vl,he + e^ v"^) of (l)-(2) on the 
interval [0,T]. We also exhibit (Propositions 1.1 and 1.2) a complete expansion as e approaches 



,1 „2^ 

,1 „,2^ ^ „,a ._ f„,al „,a2\ 



for the expression {q^^v\,v^). We find (in a sense to be specified later) that q^ ~ g" and v^ 



{vl,vl)^v'',:={vf,vf)With: 

N+l N+l 



qlit,y,e) = J2 ^'<iUt,y,o), v:{t,y,9) = Y. e'{vtit,y,o) + vi{^,y,e)). (4) 



fe=0 k=0 



These expansions reveal some time boundary layer at time t ~ (recorded at the level of the 
contribution vj, (r, •) which is exponentially decreasing with respect to the variable r) together with 
some mean evolution behaviour (described by w^). A noticeable aspect is the production of some 
dissipation when looking at the transport equation (12) on v^. The present approach is not in the 
continuation of usual k — s models [ ]. But, in the same spirit, it confirms (and justifies) that 
the interaction of oscillations can indeed be described at a macroscopic level by the introduction of 
some turbulent viscosity. 

1.2 The functional framework 

1.2.1 Sobolev spaces 

Here K denotes M, T x M or R'^. Let a = (ai,a2) € N^ The length of a is |a| := ai + a2- The 
notation 9" is for the differential operator dg^d"^. 

- Given m G N U {+00} and peWU {+00}, recall that W^-p is: 

^m,P ._ |j ^ ^P(j^) . QCf g LP{K), \a\ < m} , H" := W""'^ . 

When m e N, the space W"^'^ can be equipped with the following semi-norm and norm: 

771 

For s e M+ \ N, we can still define spaces VF^'P and iJ " by interpolation theory. 

- Let (to, n) G N^ with n < m. Define the functional spaces: 

W™'" :-{/;/£ C^([0,T];Ty'"-^'-), Vj G {0, . . . ,n}} , 
H™'" :={/;/ G C^([0,T];i/™-^), Vj G {0, . . . ,n}} , T G M+ U {+c^} , 

which can be seen as Banach spaces when provided with the norms: 

n n 

||/IU^.„:= sup E ll^t/(^'-)llw".-.-.~, ll/ll«'".":- sup E H^t /(^' Oik— 

te[o,T] ^.^0 ^ te[o,T] ^.^0 



- In order to deal with functions f(t, •) defined on M_|_ x K , which are exponentially decreasing 
in the time t G M+, and which take their values in the Sobolev space H'^, define: 

£!:^{f: sup {e''\\f{t,-)\\Hs^K))<+^}, SeRl- 

^ te[o,+oo[ -' 

- Finally, introduce £^ := f] £j, V^'° := f] V^° and V^ := f] V^^ where V € {H,W}. 

1.2.2 Families of functions 

In this paragraph, we fix some Sq G ]0, 1] and look at families of the type {fe}ee]o,ea]- 

- Assume that /^ e W"^'P{K) for all e E ]0, eo]. To control the size of f^, we can use the following 
weighted anisotropic semi-norm and norm: 

m 
(i.e) QgN2, \a\=m k=0 '^'^(i.e) 

We will say that {fe}s is bounded in WJ^'^-. when: 

III/-||Ih'"'V ^^P ll/ellw-.P < +00. 
< •■' £e]0,eo] <'■ = ' 

- Assume that /e € V™'" for all e e ]0, Eq] where V = W or V = "H. To control the size of /g, we 
can use the following norms: 

L"/2J Ln/2J 

IIMI^™.. := sup E ||£2-'a,Ve(i,.)l|y'"-, ll/ellv".;" := sup J2 \\^''^tfe{t,■)\\y^-.■ 
We will say that {fe}e is bounded in V™j" or in V™/" ^ when we have respectively: 

|||/.|||v™." :== sup ll/elly™," < +00, |||/.|||v^;? := sup ||/e||y™,^ < +oo . 

Classical embedding : for s > 1 we have H%T x R) ^ W^'^^iT x M) = L°°(T x M). When taking 
into account the dependence on e G ]0, 1], there is a loss of powers in e. Retain here that: 

3CeRt; \\fe\\L^^T^^^<Ce-'/'\\M\H^^^^^iTxR), Vee]0,£o]. (5) 

1.2.3 Decomposition of a periodic function 

Any function u E L^(T;M) can be decomposed as: 

u{9) = {u) + u* (0) , {u)eR, u*eL^{T;R), (u) = Ilu := u{9) d9 . 



In what follows, given a symbol V G {Ty"'^, iJ", W™'", H^, £|}, we will manipulate functions 
f{t,9,y) € V(T X R). We will often decompose / into its mean and oscillating parts according to 

u{t,9,y) = {u) (t,y) +u*{t,9,y), 



with, 

///(t, v) := n/(t, y) := {f{t, -, y)) G V// := nV(T x M) , (6a) 

/^ (i, 0,y):^r{t,e,y)eV^:={I-n)V{TxR). (6b) 

To signal that we consider functions f{t, 9, y) which do not depend on 6* G T (11/ — /) or whose 
mean value is zero (11/ = 0), we will use respectively (as above) the marks / and _L. By extension, 
when dealing with some operator P, we will note 

P// := PU , P^ := P{I - n) . (7) 

Be careful, in the case of operators, the composition by 11 and / — 11 is put on the right. 

The derivative de acts in the sense of distributions on the space i^(T;M). We find: 
/C :=kera0 = {u = c; ceM}, IC^ := {]<.ct de)^ = {m G L^{T-M.)- Hu = O} . 

The action dg has a (right) inverse dg^ : JC^ — > JC^ n _ff^(T;M) which is given by: 

dg\{9):^ u{s)ds- / u{s)dsde, VdeT. 
Jo Jt Jo 

1.3 Main statements 

Since we impose v ^ M > 7/2, the equations on the components u^ and u^ can be considered as 
being partially decoupled from the equation on p. Up to some extent, we can first deal with: 

£1 (e, q,, V,) := dtvl +e~^h dyv] + e^'-^ [v] dev] + e v^ dyv]) + C e^ -m-2 ^^ g^^^ „ j5i^^ ^ 
£2(6, Qe, V,) := dtv^, + £-1 h dyvl + e-2 Qgh vl + e^^-2 ^^1 Q^^2 _^ g ^2 Qy^-^ 

and then look at the remaining part as a transport equation on q^ : 

£o(e, fe, V,) := dtq^ + e'^ h dyQ, + e^'-^ {vl dgq, + e vl dyq,) + C e^'-^ q, {devl + e dyV^) . 



In what follows, the results will be expressed in terms of the quantities g^ and v^- Of course, 
the expression 

(p„ ul,ul) = is'' q,,e'' vl,h + e^" v^) (8) 

is a solution of (3) if and only if: 

/:,(£,<Z„«,) = 0, Vj €{0,1,2}. (9) 

1.3.1 Construction of approximated solutions 

We start by constructing approximated solutions for the system (9). The first step is to look at the 
two last equations of (9), where the 0(e^'^~*'^~^) <C 1 contributions (implying g^) are neglected. 
Thus, we start by considering the system: 

f /:j(e, We) := dtv\ + e'^ h dyvl + e*^-^ ^„i g^^ ^ ^ ^2 g^^i^ _ .pi^^ ^ 

\ £^(e, ve) := dtvl + £-1 h dyvl + e-^ dgh vl + s'^'^^ {vl dgv^ +evl dyV^) - V^v^ . 



Proposition 1.1. Fix an integer M G N with M > 2. Choose any integer N Cz N and any decay 
rate 6 e]0,/i[. Select any functions v^ G H°^{T x R;R'^) indexed by k e 10,N + 1}. There are 
functions 

vie f] H^, vlesr. fcelO,7V + i] 

such that the family {w"}e defined as indicated in (4) satisfies the following conditions: 
i) At the initial time t = 0, the trace w"(0, •) is prescribed in the following way: 

N+l 

<(O,0,y)=5] e^-«O(0,j;). (10) 

fc=0 

ii) For all time T G R*^, for all m G N, the family je""'^ >C"(e, w")} with j — 1 or j — 2 is 
bounded in H™^ in the sense that: 

sup sup max ||e"^/:"(e,'y^)|| < +cx) . (11) 

Furthermore, the expression nti| is determined through an equation of the form: 

djivi - (;, + i n {{d.^hf)) dyyiivi = Sk (12) 

where the source term Sk depends only on the w* with j < k — 1. 

To complete v^ into some approximated solution {q^,v^) of the complete system (3), there 
remains to identify the pressure component q^. To this end, we are satisfied to solve directly the 
transport equation Co{s, 9", u°) = where v^ is adjusted as in Proposition 1.1. Since the expression 
v^ is a function of the scales of time t and -|-, that goes for q^{t, y, 9) too. In what follows, we will 
not need to precise the way by which g° depends on the different time scales (t, -, ^, •••). 

Proposition 1.2. The context is as in Proposition 1.1. Note {v^}^ the family issued from the 
Proposition 1.1. Select functions ql € H°^{T x R;M.^) indexed by k e 10,N + 1]. There are 
functions q| with: 

{ql},e fl n^-\ fcGlO,7V + ll, eG]0,l] 

such that the expressions q'^ defined as indicated in (4) are solutions of the Cauchy problem: 

N+l 

£o(e,g?,<) = 0, ql{O,0,y) = E ^' ^"iO,y)- (13) 

fe=0 

Moreover, for all time T G M*^, for all m £ N and for all fc G |0, iV + 1], the family {ql.}e is bounded 
in T-C^'fi \ in the sense that: 

sup sup lk|(t, Oll^™ rxxR) <+oo- (14) 

£e]o,i] te[Q.T\ '-^■'^ 
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Coining back to £i and £2, we can now make the following statement. 

Proposition 1.3. Select m, M, N, i^ E N satisfying: 

M>2, m > 2 and 2v - M - 5/2 - (to + 1) - iV > 0. (15) 

Note {Wgle and {q'^}e the families obtained with Propositions 1.1 and 1.2. Then, for all j G {1,2}, 
we have: 

sup sup \\e-^ Ljie^q^^^vDW , < +00. 

eG]0,l] *e[0,T] *■ ' 

1.3.2 Existence and stability result 

The parameter e e ]0, 1] being fixed, the local in time well-posedness of the Cauchy problem (l)-(2) 
is standard, with corresponding solutions 

{Pl, <) = (e'' g|, e^ vl\h + e^' vf). (16) 

It means that, for all e G ]0, 1], there is a time T^ e M^ (eventually shrinking to zero when e goes 
to zero) such that (g|,w|) with v^ := {v^-^,v^'^) is a solution of (9) on the time interval [OjT^] with 
initial data as indicated at the level of (10) and (13). 

Fix any ReN. We can always define on the strip [OjT^], two functions q^ and v^ through the 
identity 

(g|,<) = (<Z^<)+£^(9f,z;f). (17) 

Two questions are solved below: the existence of exact solutions of (l)-(2) on a time interval [0, Tc] 
with Tc € K+ independent of e g]0, 1] and the production of controls on {q^,v^) showing that 
{q^,v^) gives indeed some good asymptotic description of (q|,i;|) on [OjTc]. 

Theorem 1.4. [Existence and stability] Assume A < 4/i. Let m, v, M, N, R E N satisfying 

M>7/2 and w,„ := min (2z/ - Af - 5/2 - (m + 3) - i?,iV - i?) > . (18) 

Let T^ the lifespan of the solution of the Cauchy problem (l)-(2). Then there exist Tc > and 
Scrit > such that 

yee]0,Sc„t], T,>Tc. 

Furthermore, the approximated solution *((7°, v^), constructed thanks to Propositions 1.1 and 1.2, 
is a relevant expansion for the exact solution ^{q'^,vf) (associated with *{p'^,u'^) threw Equation (IG)/ 
in the sense that the remainder *{q^,v^), defined in (17), satisfies the following statements. 

i) The family {qf}e is bounded in "H™ ^' -,.' 

sup sup lUf (i, •)|| ™+3 < +00. (19) 

eG]0,ee]tG[0,Te] d'^) 

ii) The families {Wg^je and {ev^^}^ are bounded in T-C^ ^ ' : 

sup sup ||"e^(i, Oll^m+a < +00, sup SUp jje ^^"'^(i, ') || //m+3 < +0O- (20) 

ee]0,ee] tG[0,Te] eSjO.ee] te[0,T<:] 



1.4 The context 

1.4.1 Historical comments 

Let N £ N* . Consider the following scalar equation of evolution: 

5t/s + ^Mv)-V,/e + ^ Q/e = 5(t,x,v), /at,a;,v)eM (21) 

where h : M.^ — 5- M.^ is some smooth function, Q is some linear operator acting on L^, and S{t, x, v) 
is some function depending on the variables (i, x, v) G ]R_|_ x K^ x M^. The unknown is the function 
f^{t,x,v). Depending on the choice of Q, the Equation (21) can be the neutron equation [ '], the 
Fokker-Planck equation [!)] or the Boltzmann transport equation [15]. In this context, it is well- 
known that the family {fe}ee]o,i] has a weak limit , say /o as e goes to 0. In general, the expression 
/o satisfies an equation implying a drift- diffusion term of the form — diva;(_D V^; • ) where D is some 
squared matrix depending on the data. The proofs of the related statements rely strongly on the 
structure of the collision operator Q which is either a bounded operator or a self-adjoint operator 
on some weighted version of L^ . 

When the operator is less regular or when there is a lack of symmetries [ ], the convergence 
concerns only the mean value p^ with respect to v, called the density. For some function qq satisfying 
adequate restrictions, we have: 



Qeit,x) := / fe{t,x,v) dv — > goit,x) . (22) 

When looking at the structure of £° :— *(£",£2)i there is some analogy with (21). Indeed, the 
expression C^^e^v^) can be decomposed into: 

C''ie,v,):^dtv, + -Tve + \ Qv, + CC{e)v, + e^'-^ NC{e,v,) (23) 



£ e 



-fidge \ ^_fhdy-\dee 

deh -fidee) ' '^ \ hd^ 



with: 



-\dey ~{fi + e\)dyyj ' \ ^ ej ■ \v^dgvj + evjdyV^ 

There are many analogies between (21) and (23). In both cases, the hierarchy with respect to 
the negative powers of e (namely e^^, e^^ and e'^) is the same, with in factor operators sharing 
analogous structures. Also, the mean value operation (ug) is when considering (23) what replaces 
the integration with respect to v at the level of (22). However, there are two important differences 
when comparing (21) and (23) : 

- The Equation (23) is a system {v^ G R^). When dealing only with the singular part e~^T-l- 
e~^Q, this problem can be circumvented by first solving the equation on v^ and then by plugging 
the result into the equation on v^. However, once the influences of the contributions CC or MC are 
incorporated, such strong decoupling is no more available. When dealing with the full system (23), 
the discussion must necessarily take into account vectorial aspects. 



- The operator Q of (23) is neither selfadjoint nor bounded (on L^). Up to some extent, it can 
be viewed as a non selfadjoint perturbation of the selfadjoint action —fj,dgg I. Still, we can compute 
the point spectrum ap{Q) of Q : L'^{T-m.'^) — > H-'^{T;R'^). We find that: 

ap{Q) := { (5 e C ; Q- S I is not injective } = {^J.n'^ ; n£N} . 

From the point of view of central variety theorems, the presence of a (point) spectral gap between 
the eigenvalue and the other (positive) eigenvalues indicates that there is a separation between 
two types of behaviours in time, a slow one and a fast decaying one, for instance in the spirit 
of [10, 1-^)]. Of course, such a separation is due to the presence of —jidgel inside Q. Again, the 
influence of this dissipation term is what relates (21) and (23). 

In other respects, singular systems like (23) have been studied in a purely hyperbolic context, 
that is when /^ = A = 0. Then, the discussion is based on tools coming from supercritical nonlinear 
geometric optics [1, 4, 6]. 

The asymptotic analysis of (23) under the assumptions retained here is clearly at the interface 
of what is done in [2, 8, 1-^)] and [1, 4, G]. Nevertheless, it needs to develop a specific approach which 
is the matter of the current contribution. In the next paragraph, we give a few indications of our 
strategy. 

1.4.2 Heuristic description 

Our analysis of £° is based on a discrete Fourier decomposition with respect to 6* e T. We can 
expand h as well as v =* (ti^,w^) into Fourier Series: 

fegZ* feGZ 

Introduce the following linear map: 

n :L2(T;M2) — ^ L2(T;M2) 

feez* 

The application 11 is clearly a projector onto the kernel of Q. Retain that: 

n o fl ^ fl , fl L^ := ker Q , dim (ker Q) = 2 . 

• a • To understand the action of the several operators in £° defined in (23), a first approach 
is to consider the simplified equation: 

dtv,+e-^Qv,^{), «aO,0) = ^%(O)e'"^ (24) 

fcGZ 

The corresponding solution u^ — *{vl,v'^) involves components v^ which can be put in the form 

feez 



For fc = 0, we find that VQ{t,T) = WqIO) ^iid: 

The second (constant) term in VQ{t, r) is in general non zero and it conies from contributions inside 
^^(O, •) which are polarized according to (/ — 11) L^. Thus, even if v^{0, ■) presses only on the positive 
point spectrum, the corresponding solution Ve is not necessarily exponentially decreasing in time. 
We can see here a first effect of the nonself adjoint part inside Q. 

For keZ*, noting H := {p e Z*;p ^ k,p ^ 2k}, we have vl{t,T) = ^^(O) e^'^'^''' and: 
vl{t,T) = -ihkfi-^k-^i^{0)^2ikh2k wifc(O) r 6-'^'="^ 

By bringing together all constant terms (in t) inside an expression v^{t,9) which here does not 
depend on t, these formulas fit with a decomposition like (4). 

• b • Next, consider the more elaborated model: 

dtVe+e-^Tv,+e-^Qv,^0, i;,(O,y,0) = ^ 1)^(0,2/) e^'=^ (25) 

feez 

One can expect that the intermediate singular term e~^ T produces the scaling t/e. However, 
such an effect does not appear here. On the one hand, the contributions polarized according to 
(/ — n) L^ are mainly handled as in paragraph a. On the other hand, the H L^ parts disappear by 
a combination of two arguments: 

- Due to the relation J^ hh'd9 = 0, we can use the following algebraic identity: 

noron = o. (26) 

- We can absorb the extra term (/ — n)Tn through some ellipticity inside Q. Indeed, in what 

follows, we seek v^ as an expansion of the form v^ — wo + ewi + O(e^). Assuming that 
Wo = Huq G ker Q, we can observe that: 

{I^U){e-^T + e-^Q){vQ+evi) = e'^ [(/ - fl)rflt'o + (/ - ff)Q(/ - 5) ^i] +0(1). 
Now, the idea is to adjust vi conveniently in order to remove the 0{e^^) contribution. 

In practice, the implementation of these arguments must be done with care because the different 
terms which come into play are more tangled than what is indicated above. 

Note that a normal form approach (in the spirit of [.,]: meaning to change v into (/ + eM)v for 
some well adjusted operator M), can be tried to get rid of T. However, such a method seems not 
to succeed. There are always remaining 0{e~^) terms and, all things considered, to deal with the 
actual diagonal form of T appears to be more suitable. 
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• c • Finally, consider the full system (23). Our aim is to describe the asymptotic behaviour of 
the family {weje on a time scale of the order i ~ 1. To this end, we have to understand the 0(1) 
contributions brought by the singularity e~^ T + e~^ Q. This singular term is a perturbation of the 
self adjoint operator Qo '■= £~^^J■^eeI■ This perturbation is of two types. 

- The interactions (at order one) between Qq and e^^T turns out to be the source of some cre- 

ation of diffusion. The mechanism is similar to the one met in the drift- diffusion phenomena. 
Moreover, Qo + e~^T is a diagonal operator. The components of the velocity are decoupled 
and the discussion deals more with scalar arguments than with vectorial arguments. 

- We perturb Qo + £^^T by I r, , n ) ^^ order 0. A first effect is that Q is not selfadjoint. It 

also induced some strong coupling at order between the two components of the velocity. One 
aspect of the construction is to prove that this strong coupling do not disrupt the production 
of dissipation. The discussion has to take into account vectorial aspects and one issue is to 
match the initial data between the slow profile w| and the fast profile v^ . 

Moreover, we have to determine the effects of CC and J\fC which are of two types. First, the 
presence of CC and AfC reinforces the coupling. Secondly, it induces nonlinear interactions which 
are delicate to deal with. In particular, in the critical case M = 2, it becomes necessary to exhibit 
transparency phenomena in order to achieve the analysis. 

In this article, we propose (Proposition 1.3) and we justify (Theorem 1.4) a complete expansion 
for the family {feje. It is the occasion to analyze precisely the linear features and the non linear 
aspects alluded above. 

1.4.3 Heuristic arguments for the energy estimates 

The variable *(g^,t;^) can be interpreted as the solution of the linearized operator C at *(g°,w°) 
perturbed by some small non- linear terms (when i/ and M are large enough). It can still be 
interpreted as the solution of the linearized equation of System (3) at point *(0, 0, h) + {e"^ q^ , e^ v^) 
(perturbed b some small non-linear terms). To underline the difficulties to obtain estimates on a 
strip independent of e we only consider estimates over (pf , ul) solution the linearized system (3) at 
point *(0,0,/i(6')): 

%pi _ = So, 

'hdyuf -VJ4 =5i, (27) 




hd„u 
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2f,u ,,U 



e dgh u 



'2: 



for some sources S := '^{So,Si, S2) in H°°{T x M ; R^). It is a parabolic-hyperbolic system singular 
in s. 



Purely hyperbolic approach. We first consider that /i = A = so that the dissipation vanishes. 
We perform classical L^-estimates on Equation (27). We obtain: 

\\{pl,ui){t,-)\\^,<e''^' sup \\S\\^,, with C,<C{l+e-^\\dgh\\L^). 

ie[0,Tj] 
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Yet, it indicates that classical energy estimates only provide a control over the solution for time of 
order e^ . In particular for bounded time the solution can exponentially increase with the time t. 
Furthermore, let us consider the singular transport equation: 



dtv + eh{0)dyv ^ 0, 



^=0 = ^'0 ■ 



The solution is explicit v{t, 0, y) = Vq{9, y — e^^th{6)). In particular each time we differentiate with 
respect to 9, we loss a power of e. Thus, classical Sobolev space are not well suited for the control 
of this family of solutions. We have to introduce anisotropic Sobolev spaces (defined page 4) for 
both the velocity and the pressure. 



Parabolic-hyperbolic approach. To go further in time, we have to consider the dissipation. 
The operator Ve, is positive and satisfies some coercive estimates. There exists a positive constant 
c such that for any function / € H^{T x K), 



V£e]0,l], -{Vefj)>c[\\e-'def\\l,^^^, 
It has two consequences. 



\dyf\\ 



L2(T> 



*.(v,/). 



(28) 



- At fixed £, we should obtain a regularization of the solution. The velocity ul is in L'^Hg (see 

Inequality (101)). This is the regularization phenomena. 

- Considering the dependency in e, the dissipation should also absorb some singular terms. First 

the squared term can be estimated as follows: 



£ ^ / dghvl^fdOdy 

TxR 



< 



m 



|,'||2 



\e-'devi\\ 



ellL^ 



Thus it only seems to be singular of order one (in e) instead of being singular of order two 
(in e). Furthermore it is also desingularizes the singular transport. We can obtain estimates 
over the velocity in the classical Sobolev spaces whereas the pressure is still estimated in the 
anisotropic Sobolev spaces. It indicates that the estimates over the velocity and the pressure 
have to be done separately. 

Here, the addition of the parabolic aspect in the discussion still does not allow us to obtain a 
control over (p^,u^) for time of order one. Some additional arguments are required. 



Singular change of unknowns. To keep on desingularizing the term e ^dghul^ we consider 
the change of unknowns: 






ql, K := wf , 



~:if. 



,2£ 



The system (27) becomes: 

dtuf 



+e~^hdypl 



u 









So , 
eS2. 



(29) 



where the operator Q^ is defined in Equation (97). We can notice that: 
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- The term e ^^e/iu^^ is desingularized into e ^dghuY- 

- However, the dissipation is turns into the operator Q^. It can no longer satisfies Inequahty 

(28). Assuming /i is large enough, it is still true (c./. Lemma 3.5). 

Thus performing the same estimates for system (29) as the one done in the previous case should 
lead to a control over (gf, ul) in L^-norm for time of order one (t « 1). 



Conclusion. In Section 3, we justify that those heuristic arguments work for the complete Sys- 
tem (94). Some technical arguments must be added to deal with the complete System (94). Indeed, 
it is obviously nonlinear and the pressure and the velocity are coupled. Nonlinear terms has to be 
studied carefully. 

As indicated, the pressure is expected to be controlled in anisotropic Sobolev spaces whereas 
the velocity is estimated in the classical Sobolev spaces. Thus it indicates that we have to deal 
with the problem of the velocity and the pressure separately. However those variables are coupled 
by the term 

-^—^ \de, edy) {q^ + e\ff , 

in Equation (94). The constant v has to be large enough so that the pressure does not interfere 
too much with the velocity. Of course an other issue is that the pressure is only estimated in the 
anisotropic Sobolev spaces. We can go back to the classical Sobolev spaces using the equivalence 
of norms (92). It has a cost in power of e for each derivatives to estimate. It explains why we lose 
(to + 3) precision in the definition of Wm (see Equation (18)). 

1.4.4 Contents 

What follows is divided in two main parts: Section 2 and Section 3. 

The Section 2 is devoted to the construction of the approximated solutions {q^,v^). The first 
step is to show the Proposition 1.1. 

- In this purpose, the paragraph 2.1 deal with the velocity field u°, that is with the equation 
£°'{e,v^) = 0{e^). The limit case M = 2 is special because in this situation the non linear 
terms can interfere at leading order. 

- In the paragraph 2.1.5, we are able to exhibit the control (11). 

The pressure component q^ is incorporated at the level of subsection 2.2. Then, the complete 
construction of {q^, v^) can be achieved in the form of Proposition 1.3. 

The Section 3 is concerned with energy estimates. We first state the Proposition of control 
of the velocity and the pressure. In particular we deduce estimates stated in the Theorem 1.4. 
In the subsection 3.1, we look at the equations Ci and £2- To this end, we crucially need the 
properties brought by the dissipation. In the Subsection 3.2, we inject the informations which have 
been obtained at the level of Cq- By this way, we can deduce controls concerning the pressure 
component. 
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2 Construction of the approximated solutions 

We recall that M is assume to be larger than 2. This section is dedicated to the proof of Propositions 
1.1, 1.2 and 1.3. 

2.1 Approximated velocity 

In this Section, we construct an approximated velocity. Since M > 2, it follows that non linear 
effects are present. We are forced to work with the two time scales t and e~^ t together. We 
construct expansions, 

N+l N+1 

fc=0 fe=0 

and plug the expression -y| + v^ into £" at the level of (4). We obtain: 

dtvUt, .) + e-'hdyvlit, .) + e*^-2 {vl^t, .)devlit, .) + sv'/{t, .)dyvl{t, .)) 

+ e-^ ' {0,dghvl%t,.)) -V,vl{t,.) 
+ £-2 drvf{t/e^, .) + e-^hdyvf{t/e^, .) + e*^-^ yif(^t/e^, .) devf{t/e^, .) 

+ e''-'v^Jit/e^,.)dyv{it/e^,.)+ e-^ ' {0,dghvlf{t/e^,.)) -r^vlit/e^,.) 

+ e*'^-2 {v^t, .)devi{tle\ .) + evf{t, .)dyvi{t/e\ .)) 

+ e^'-^vlf{tle\.)devl{t,.)+e''-'vlf(tle^.)dyvl{t,.). (31) 

Fix any time T E M.*^. Define 

C'^'ie, vl) := dtvl + e-^hdyvl + e*^^^ ^^I^q^^s^ ^ evfdyvl) + e'"' * (O, dehv]') - V,vl . (32) 

Recall that v^ £ £g° is assumed to be exponentially decreasing with respect to r e M+. Since 

^-st/e' = 0((£Vi)°°) = 0{e'^) , V (t, N) e ]0, T] X N , 

when looking at the Equation (31) for times t g]0,T] with in view a precision of the size 0{e^), 
all terms involving v| can be neglected. Now, the idea is simply to extend this (relaxed) smallness 
requirement on the whole interval [0,r]. Briefly, we seek w| so that 

C-%e,vl) = 0(6^), te[0,T]. (33) 

The Equation (33) can be completed with some initial data 

N+l 

<(0, 0, y) = vf{e, y)^J2'' ^fc°(^' y) ■ (34) 

fe=0 

Clearly, it suffices to specify vf to determine what is v^{t, •) for t £ [0,T], by solving the Cauchy 
problem (33)-(34). Now, in order to select vf conveniently, we have to take into account what 
happens for small times, in a boundary layer of size e^ near t = 0. To understand why, just come 
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back to the study of (31) for i ~ e^ or t ~ 1. Then, the contributions brought by u| can no more 
be neglected. Considering (31) with the information (33) in mind, it seems natural to impose 

C^f{e,vi,vl) = 0{e''), r e [0, 1] (35) 

where we have introduced 

£"/(£, viX) := e-^drvi{r, .) + e-^hdyvi{T, .) - V.vHr, .) 

+ e''-' (vlHr, .)deviiT, .) + ev'Jir, .)dyvf{T, .)) +* (O, dehvlHr, .)) 
+ s''-' {vl^{s\ .)devi{T, .) + e«f (eV, .)a,«/(r, .)) 



- £ 



Af-2 



(vlf{r, .)dovtie\ .) + ev'Jir, .)dyvlie\ .)) ■ (36) 



Assume that the data z;|° is known. The Cauchy problem (33)-(34) furnishes w|(i, •) for t S [0,T]. 
In particular, it gives access to all derivatives (9t)'i'£*(0, .) with I £ N. Therefore, we can go further 
in the analysis by replacing in £°'^(e, i;^,t)|) all the expressions w"(e^T, .) by their corresponding 
Taylor expansions (say up to the order A^ — 1) near t = 0. As long as t e M is fixed, this operation 
is justified. From now on, we look at 

£'^^*(£,w|') = 0(£^), TeR+ (37) 

where the definition of C^-^* is 

£"^*(e, vl) : = £-2 g^^^(^^ ) ^ e-^hdyV,{T, .) - P,v,{t, .) 

+ e''-' {vlir, .)dev,{T, .) + evl{T, .)9,z;,(r, .)) + s'^ * (O, dghvl{T, .)) 




{dtYvl^iO, .)dgV,{T, .)+eJ2 '-^{dtYv'/iO, .)dyV,{T, .) 



1=0 



N-1 ,^2^\l 



^^{dtYdevl{Q,.) + evl{T,.)Y,^^{dtYdyvl[0,.)\ . (38) 

1=0 ' 1=0 ' J 

To be coherent with (4), we have to impose 

«/(0,.) =«f(.) := (v^-vtM-), (39) 

where w"(0, •) and wf (0, •) are prescribed as indicated in lines (10) and (34). Now, we can recover 
some u/(t, •) for t G M+ by solving the Cauchy problem (37)-(39). 

The difficulty comes from the condition v^ G £g° . Nothing guarantees that the criterion v^ G £g" 
can be verified for some well-chosen u|". To show the existence and the uniqueness of such a data u|° 
is in fact what matters. The extraction of an adequate function vf' is clarified in the construction 
described below. 

For the sake of brevity, for k > M, introduce the following notations: 

J{M,k):^{{i,j)elO,N+lf ; i + j=k~M} , (40) 

X(M,A;):={(i,j,Oe[0,iV + lf X 10,N - Ij : i + j +21 = k - M} . (41) 
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For fc < M, we set J{M, fc) = and Z(Af , k) = 0. We also adopt the conventions v^ = vl = for 
k — —3, k ^ —2 and k = —1. Let us now go into the details of the BKW calculus. 

The first step is to inject some expansion w| like in (30) into the Equation (33). By this way, 
we can obtain a cascade of equations concerning the unknowns vf. = *(wj ^£,'^2 fc)- More precisely, 
for k G |0, iV + 3], we have to consider 

= ^J■^eevtk + fJ'9yyvt,k~2 + ^^eew^_i + A^ej^w^ j.^2 . (42a) 

{i,3)eJ(M,k) (tj)eJ{M+l,k) 

= IJ-degvlk + fJ'dyyvlk-2 + ^doyvlj^^^ + \dyyvli^_j . (42b) 

The next step is to plug some expansion v^ like in (30) into the Equation (37). By this way, we 
can obtain a cascade of equations concerning the unknowns vl — *(^i ^1^2 fe)- More precisely, for 
fee |G, A^ + 1], we have to consider 

{i,j)eJ{M,k) (ij)ej(M+i,fe) 

+ E {<'^ (^*)' (devi^m + id,y (vumdovQ ^ 

{i.,j,l)GX(M.k) 
{i,jU)eX(M+l,k) 

= /^^eei'f.fe + f^dyyv(,k-2 + ^9eev(^k-i + •^^02/«2,fc-2 ' (43a) 

{t.'j)eJ[M,k) {i,j)eJ(M,k+i) 



r' 



+ E {4Adtf {devl^m + {dt)\vl,mdevQ 

(ij,OGl(M,fe) 

+ E [4Adtf{dyvi^m + {dt)\vi,mdyvQ 

{i,jU)eX(M+l,k) 

= ^^^eevlf, + f^dyyv^j^^^ + Xdeyv(,^_^ + Xdyyv[j^^^ . (43b) 

In view of (10), we can associate (42) and (43) with initial data vf and v^ satisfying the restriction: 

vl{0,e,y)+vl{0,d,y) =vt{d,y), V/c G lO,iV + 1] . (44) 

Proposition 2.1. [Solving (42) and (43) together with (44) o-nd the condition v^ G £g°] Fix a time 
T e M;, a number S € ]0, fi[ and, for all k € [0, ..., N + Ij, functions v^ € H°^(T xR). Then, the 

16 



conditions (42*^), (43i^) and (44) have a unique solution such that 

(vl,vi)en^xEr, VfcelO,iV + ll. (45) 

Moreover, the component v^ can be identified through the homogenized Equation (12). 

The proof relies on some induction on the size of N, based on the following hypothesis of 
induction: 

'HJ^{N) : " The Proposition 2.1 is verified up to the integer N". (46) 

To go from A^ up to iV + 1, we will need a succession of lemmas which are produced below. 
Before going into the details of the analysis, we give below a brief description of what happens 
depending on the choice of M. 

- When M > 3; the non linear terms are rather small, the construction is somehow linear. 

- When M — 2, the non linearity becomes critical and a few arguments must be added. For 
instance, if we write the Equation (42a) for fc = 2, we can notice a Burgers' term 

5tUi + f^ 9yVii + vIq devfg = ^ deevl2 + M dyyvf^ + Xdeevl-^ + A^sywlo (47) 

and also two contributions dggvf ^ and dggvf ^ to be calculated in function of vf q, with ap- 
parently a non linear dependence with respect to vf q. 

- When M — 2 again, another effect of the non linear interactions is the apparition at the level 
of (50) below together with (62), of a source term S^ which can depend on v^. 

- However, there are transparency phenomena at work which come from the initialization pro- 
cedure. Indeed, knowing that v^ = for k e {—3, —2,-1}, the Equation (42) in the case 
A; = and M >2 reduces to fidggVQ = 0. In other words, we have to impose 



,s-L 



(/-nXo^O. (48) 



It follows that vf Qdgvf Q = 0. In the same way, all other apparent non linear contributions 
will disappear. Therefore, the remaining term Ilvf q can be determined apart without seeing 
any non linear effect. 

2.1.1 Technical Lemmas 

A consequence of the (point) spectral gap. The system (43) is made of two evolution equa- 
tions of parabolic type, based on dr — i-J^dgg. This falls under the following framework. 

Lemma 2.2. [Fast decreasing under a polarization condition/ Let m G N and 5 G]0,/i[. Select 
Wo e iJ™+2(T X M) and Sq G (£"+^)-L(T x R), that is such that USq = 0. Consider the initial 
value problem: 

drw - ^ deew = Sq , w\^^g = w^ . (49) 

For all T G M+, there is a unique solution w G C"'^([0, T]; H"^{T x M)) to the Cauchy problem (49). 
Moreover, if the initial data is well prepared in the sense that H-Wq ~ 0, then IIw ~ for all t G [0, T] 
and wGff (T x M). 

The proof of Lemma 2.2 is very easy. It will not be detailed. 
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ql// 


I qnl/l 


+ Pi^{I- 


-n> 


i^b. 


ql// 


+ </ 


+ Pt{I- 


-n> 


^b 




+ Py/n 


ylk + ^2^ 


U- 


n)«i 



Interpretation of the system (42). The Lemma below is intended to look at the system (42) 
otherwise. Indeed, there is a difficulty when dealing with (42) since the knowledge of v1_2 seems to 
require the determination of vf._i and v^, that is the identification of terms w^ with indices j greater 
than fc — 2. An important remark is that such a dependence can disappear when the projector 11 
is applied. The possible dependence in the terms u^ are recorded in the source terms S^ (with 
explicit formulas). For this occasion, the two cases M > 3 and M — 2 must be distinguished. This 
fact can be formulated in the following way. 

Lemma 2.3. [Non-linear homogenizationj Assume that the functions v| with fc G |0, A^ + 3] are 
solutions of the system (42). Then, for all fc G JO, A + 1], the part Tlv^ = *(nuj^,nw2fc) is a 
solution of 

mvi, - (m + ^ mOo'hA dyyUvi, = S!^/, + S-\l + Pi^(/ - JX)vl^ (50) 



In the above system (50), the four operators P^^, ^2^; P2 '^^'^ Q'2' ' '^•^ "^^^^ '^^ ^^^ source term 
Sj = *(*S'/j,, 5*2 ^.), are defined along lines (51)- ■ ■ ■ -(58). On the other hand, the contribution 
gnl// ^ t(^s^l//^snl//^ ^^ g^^^^ f^y (59)-(60). 

Proof of Lemma 2.3. The matter is to identify the contributions brought by the non linear terms. 
Let us project (42) according to V^/. Since J{M + 1, fc + 2) = J{M, k + 1), this yields 

dtUvi, + n {hdyvi,+,) + n I ^ v^, devi^ 

\(i,f)eJ{M,k+2) 

\{i,j)€j{M,k+l) J 

dtUvi, + n (hdyvi,^,) + n I Y. «M 9evl, 

\it.])ej{M,k+2) 
+ n J2 ^hdyvlA+Tl{dehvl,^^2)=fidyyUvlf^ + XdyyUvl,^_^. 

\{i.j)eJ{MM+l) J 

In what follows, we will use the system (42) and many integrations by parts in order to interpret 
n(/i dyvf j._(.i), n ( h dyV2 k+i ) and n ( Ogk vf f.^2 ) • The goal is to show that these quantities can 
be expressed in terms of the v^ with j G |0, fc]. 
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o Study of Il{hdyvf /.^i). Exploiting (42a) with k + 1 in place of k, we can deduce 



-nld,,\h)dyi Yl <,9evl,+ Yl 

^ \ \{ij)eJ{M.k+i) ii,MJ{- 



vL dyVlj 
J(M,k) 



+ \ n{hd,,\h) dyyvi^) - - n(a-/(/i) a,ee«UO ■ 

Recalling (6) and (7) , we have to deal with 

ii{hdyvi^^,) = - - n {{d.'hf) dyy^vik ~ s'l^ - p[^{i - n)vi^ 



^ V \{i,j)eJ(M,k+i) (t.j)eJ{M,k) 



with the conventions, 



S{^k ■= -- deg{h) dy (9iWi,fe_i - f^dyyvlk_^ - A9eyi'lfe_i) , (51) 

Plf :=-Il {dg,\h) deeyf) - ^ H {hd^,\h) Oyyf) . (52) 

Study of Il{hdyV2 fe+i)- Exploiting (42b) with fc + 1 in place of k, we can derive 

n {hdyvi,+,) = n(a,-/(/i) dyde9vi,+,) , 

^ V \{i,])eJ{M.k+i) ii.j)eJ{M,k) 
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We come back to the Equation (42a) in order to change the last term in this sum. This yields 
= — n {dshdgg{h) dydgg{I - n) {dtvlk-i - ^^dyy^lk-i " A^eyt-a./c-i)) 

+ 4 n {dehdg,\h) dyd,,\i - n) (hdyvi^ - xdegvi^)) 

V \{t..3)eJ{M,k+i) {h3)eJ{- 



M,fe) 



This amounts to the same thing as 



1 



^(hdyvi,^,) = -- n {(d.'hf) dyynvl, ^ s!^[ - p^'^ (/ - n)t;i% - p^'^^ nvl 



f^ 



k 



-Qi'^{I-U)vl, + ^uid^g\h)dyi Y. <.9gvl^+ E ^h9y<3 

^ V \{t,3)eJ{M.k+i) ii,])eJ{M,k) 

]^uidgh a,-/ (h) dyd^g' (/ - n) J2 ^M dovlj + E ^2,. a,«j,, 

V \ii,j)ej{M,k+i) {i,j)eJiM.k) 



with the notations 






2 dehdjih) dydgg^{I - H) (ati;^.„i - ndyyvl,^_^ - A^e^Wa.fc-i) . (53) 



^2 V := - n ( ^ a^/i a,-/ (/x) dyd,,' (/ - n) (/^a,/ - xdgef) ] , (54) 



Q'if--=--U{hd^g'{h)dyyf). (55) 

A* 

o There remains to compute Il{dgh vf 1^,2) ■ This is again (42a) with this time fc + 2 in place 
of Ac. 

= ^ n {d^'ih) {dtvi, - t^dyyvi, - xdgyii - u)vi,)) 

A* 
+ - n{dg\h) {hdyvl,^, - Xdeevl,+,)) 

'^ \ \(i,j)€j(M,k+2) (i,j)gJ'(M,fc+l) 
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We want to remove the presence of vf j.^-^. The relation II{hdg ^h) = allows to write 

Il{d^\h) ihdyvl,+, - Xdeevl,+,)) = n{d^\h) [hdy ~ Xdee) d,,^ {I - n) deevl^^,)) . 
The part dgevl j.^^^ can be extracted from the Equation (42a) with fc + 1 in place of k. We find 

+ 7:2 n {dg\h) {hdy - \dee)dg^{I - n) {dtvlk-1 - l^dyyvlk_^ - Xdeyvlk-i)) 
+ ^ n {d,\h) (hdy - \dee)d^e\l - n) {hdyvl^ ~ Xdeevlk)) 

V \{i,j)eJiALk+i) (%,3)eJ(M.k) 

^ \ \{i,j)eJ{M,k+2) (jj)gj(M,fc+l) 

This yields to 

n (dehvi,^,) ^ -5^^f - p^2± (/ - n)«j,, - pf // n«j,, - g^ (/ - n)vi, 

\ \(t,3)€JiM.k+l) 



-J ' 

{i,j)eJ(M,k) 



'^ \ \iiJ)eJ{ALk+2) 



■3 

{i,j)eJ{M,k+i) 



E "2,» ^y^i,: 



with 



Si% :=-i dg\h) {-hdy + \dee)d,^{I - n)(a,«tfe-i " m5,,«U_i - Xdoyvl^.^,)) , (56) 



A* 



P^V := J_ n (9-i(/i) {-hdy + A9ee) 9^7 (/ - n) {hdyf - Xdeef)) 
A* 

-^n{^^\h){^tf-^i^yyf)) , (57) 

Q^2/:=^n(a,-i(/i)(5e,/)) . (58) 

A* 

Briefly, for all fc e N, we denote 

"^2,^: •— ^2,k + ^2,k + -^ ^!/anW2^fe_i , Pj := Pj + ^2 J Q2 •= Q2 + Q2 • 



21 



Conclusion. Combining all informations together, we finally obtain (50) with 

^ V \(2,j)eJ{M,k+i) {ij)ej(M,k) 

-n| E <^9evl^+ E <^9y<X (59) 



and 



^ \ \{i,j)eJ(M,k+l) (r,])eJ{M,k) 

-X^uideh d^^{h) dydg,\l - n) Yl <^ 9evl, + E <^ 9y<. 

'^ V \(»j)eJ(M,fc+i) (i,j)€J{M,k) 

- -^ n [ 9,-i(/i) (hdy ~ \dee)d^^{I - n) ^ vl, devl, + J] z;^ dyvl 

^ V \(*j)e-7(Af,fc+i) (ij)ej(M.fc) 



^ V \(i,])eJ{ALk+2) 



3 

(i,j)eJ{MM-+i) 



E "2,^ dyVl 



-n ^ «M5e<,+ E ^hdyyh]- (60) 

\(»j)6J'(M,/c+2) (ij)gJ-(M,fc+l) y 

D 
The definition of Sl^" involves the sets J{M, fc), J{M, k + l) and J{M, fc + 2). Note that 

[M>3, fee {fc,/c + l,fc + 2}, (i,j) e J^(Af, fc)] =4> i < fc- 1 and j < /fc-l. (61) 

Thus, the presence of vf. inside S^ is not allowed as long as M > 3. On the contrary, when M — 2, 
it becomes effective. This remark can be formalized through the following statement. 

Lemma 2.4. [Refined description of the source term of (5U)j 

i) The expressions S^ only depend on the u| with j g |0, fc — 1]. More precisely, we have 

where f^ := *(/i ^7/2^) ^.^^ homogeneous linear functions of their arguments. 

ii) When M > 3, the expressions S^ only depend on the u* with j G JO, fc— 1]. More precisely, 
we can retain that 

qnlff _ rqnlff , s s \ 

^k — Jk v"0' ■ ■ ■ T^k-l) 

where /^" := *(/]'fe 1/2^ ) '^^^ quadratic functions of their arguments. 
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in) When M = 2, the expressions S^^ only depend on the v^ with j G |0, fc] . The influence of 
vf, can be specified through a decomposition of the form 



Q'nl// _ finW/ /„s 



rqnl// 
Jk 



, pqnljl rqnl II-. 
\Jl,k '•/2,fc )■ 



Again /^" is a quadratic function of its arguments whereas SP{vq) := *(5Pi(fg), 5^2(^0)) ^^ ^^^ 
linear differential operator defined according to 

SP,{v^o) Vk ■■= - ^(^'l^o dgvl, + vl, dev^,) - ^ H {d^\h) {v^^ dev^ + t'l,.- dgvl^)) . ^^^^ 

From now on, we simply note /^ := ff! + f'^ . 

Proof of Lemma 2.4- As already mentioned, the statement i) is a direct consequence of (61). The 
linear aspect of /j. is a consequence of the formulas obtained in (53) and (56). On the other hand, 
the quadratic aspect of /j!" is obvious in view of (59) and (60). It proves ii). 

There remains to consider the situation Hi) where M = 2. Note that J^{2, k + 2) = J'(0, k). In 
view of (61), we have to concentrate on the contribution of 

/ 



n 



E 



\ 



y{iJ)eJ{OM 



n 



M '5e«i,j 



V 



fj- 



"2,3 



\i:i,j)eJ{o,k) J J \ii,fKJ(o,k) J J 

In the sums above, only the extremal indices {i,j) — (fc, 0) and {i,j) ~ (0, k) give a contribution to 
include in SP{vq), leading to (62). D 



2.1.2 Analysis of the system (50) 

Projection of system (50). The system (50) is clearly an evolution equation of parabolic type. 
As such, it can be completed by initial data. But to solve it, we also need to identify the extra 
terms (/ — n)w| j, with • e {1,2}. To this end, it suffices again to exploit (42). Recall that 
V G {W™''P,H^,yV^rP''' jTi^jSg} and de&ie the linear continuous isomorphism 

$ : V(T X M) X V(T X M) — > V*°* := V-L(T x K) x V//(M) x V-L(T x M) x V//(M) 

(/i,/2) ^ ((/-n)/i,n/i,(/-n)/2,n/2). 

By construction, the expression Vjf :— *^vf. must be solution of the system 



AV,! 



/ t^dee 








-Pi^ 


Py 





-Ts 


-Ts 


lidee 


\-H^ 


pi// 


-Q\^ 






Pyl 



Vk 



pnlll 
Jl.k 
fnlA^ 
J2,k 
\ f^l// I 
\j2,k I 



-■ Jk ' 



V,' 



»// 

"2.k 



f-alll 



where the term /j, is given by Lemma 2.4 whereas Tg and Py are the operators defined by 
T,/ := (/ - n) {dehf) , Pyf := dtf - (a^ + -H {{dg'hf)) dyyf. 



(63) 
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Description of fj}'"^. As for /^* we decompose /^'^ into /^'^ '■— fl^ + /^" where the hnear 
part fl^ is defined by 

fuc ■={! - n) (9i«i,fe_2 - fJ'dyyvlk-2 - >^deyV?,j,_^ + hdyvlj,_^ - A9eeWi,fe_i) , (64) 

ft^k ■={! - n)(9ti'lfe-2 - l^dyyvlf,_^ - Xdeyvlk_2 + hdyvlk^i - Xdyyvli,^^) , (65) 

whereas the quadratic part f^" = * (/f^ , /|'^ ) is given by 

fJ^:=(I-U)l Y. <^^^<.+ E 4.^.<.)' Pe{l,2}. (66) 

\{i.j)eJiM.k) {t,j)eJ{M+i,k) J 

Lemma 2.5. [Description of fj}'"^] The functions f\^^ and /j?" only depend on the v'j with 
j e |0,fc — 1]. They are respectively homogeneous linear and homogeneous quadratic functions of 
their arguments (uq, . . . , v'^_^). 

Proof of Lemma 2.5. Just look at (64)-(65)-(66) together with (61). D 

Analysis of the system (63). To study the system (63), we have to distinguish the general 
case Af > 3 from the critical case M = 2. 

o When M > 3, According to the Lemmas 2.3 and 2.5, the above right hand term /^' can be 
viewed as a source term. 

o When M — 2, f^ '' is no longer a source term. Now, recall that we can exploit the condi- 
tion (48). This information is essential. It induces many simplifications when computing SP{vq). 
We find that SPi{vq) = whereas SP2{vq) can be reduced to the following linear (non differential) 
operator 

SP^iv^o) vl-~ ^«k devlo) + -, <i n (hvlX) ■ (67) 

At first sight, the expression SP2{vq) Vq depends in a non linear way on Vq. However, we can again 
exploit the condition (48) (which says that vf-{^ = 0) and then apply (67) with fc = in order to 
obtain further cancellations. There remains 

SP{v^,)v^,^0. (68) 

From now on, since there is no more ambiguity, we can omit to signal that SP{vq) depends on 
v^, and in fact only on vl^^. We will most often note SP{v^) = SP{v'J) = SP = *(S'Pi, S'Pj). 
Since SP2 ^ 0, the formulation (63) must be changed. This time, we have to deal with 

AV,^ = \K't.!tLnfjtl) . n := '-^ vl (69) 

where 



A 



1 


l^dee 








o\ 




-P[^ 


Py 










-Ts 


-Ts 


l^dee 





\- 


-Pi^ - SP2 


_pi± 
^2 


'Q\^ 


Py) 



(70) 
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The similarities between A and A are obvious. These two matrix valued operators have both a 
triangular structure. The difference, when passing from A to A, concerns only the perturbation in 
the bottom-left position (4, 1). This particularity plays a crucial part in the discussion below. To 
simplify, we present below the result in a smooth setting. 

Lemma 2.6. [Solving the system (63) or (69)/ We assume that the condition (48) is verified. Select 
a function V^f = *(Vo^^^ V^}^^) £ H°°{Rf and a source term F = *{F^-^ , F^/^ , F^-^ , F^/^) G (H??)*"*. 
Then, for all T G M^, the problem 

where V G {yl, ^}, has a unique solution V in (Hf?)*"*. 

Proof of Lemma 2.6. To solve (71), the strategy is to argue line after line. We start to solve the 
problem for the operator A. 

- First line. Since the operator ^dee '■ {T'L't) — ^ (^t") i^ invertible, we can define without 
ambiguity (and with no choice) 

V^^ := fi-^dgg^F^^ G (n^)^ . (72) 

- Second line. Observe that P^-*^ : "H^ — > (^7?) • The next component V^^^ must be a solution 
of the heat equation (in t and y) 

{PyV^" =F^'/ +P[^V^^ e{H^)" , {V^"X=o^V^" e{H^)" . (73) 

Obviously, there is a unique solution on [0, T] of this initial value problem. It does not depend on 
e. In other words, it is such that V{t, •) G {W^)'/ for all t G [0,T]. 

- Third line. Since Tg : T-U^ — > (^t") the component V"^^ can be obtained through the 
formula 

V^ = M"ia-/(F2^ + T,v^^ + T,v^/')) G [n^t ■ (74) 

- Fourth line. We can use the same argument as in the second line. It suffices to check that by 
definition P^-.U^ — > {U^)^ and Q\:n^ — > {U^)^ . Then, there remains to solve 

{ PyV^// = p2// + pjiyli + pV/^1// ^ gi±y2± ^ (^op)// ^ {V^'/X^^ = V^" G (H??)// . 

Note that the triangular structure of A is crucial in this procedure. 

To solve the system (69) that is replacing the operator A by A, the only change in the above 
proof is at the level of the fourth line where the supplementary source term SP2V^^ must be 
incorporated. D 

2.1.3 Analysis of the system (43) 

Projection of the system (43). In this paragraph, we consider the parabolic system (43) which 
can be associated with some smooth initial data w{(0, •) G H°°(T x K). Classical statements (see 
for instance [11]) say that the corresponding Cauchy problem has a unique global solution u;J such 
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that vl{t,-) e 71(1 X M) for all t e M+. The difficulty is the following. The variable r is aimed to be 
replaced by e^^ t with i fixed and e -> and, since the original equation (3) contains nonlinearities, 
we cannot allow any (uncontrolled) growth with respect to t. To get round this problem, we will 
instead require a rapid decay when r — >■ +(X) but this necessitates vl (0, •) to be selected conveniently. 

To see how to adjust ^{(0, •), we can interpret (43) in the form 



BV^:^ 



fPe 0\ /<A 

nl// 

9i.k 

nl± 



dr 

Ts T, Pg 
\Tf drj 



.92X 



V, 



f ttf,(nJ\ — 



*K) 



/J- 



(75) 



where Tf and Pe are the operators defined by 



Tff:^nidghf), Pel-.^drl-tideel. 



Description of g^'^. We decompose 5^?'^ into is linear and quadratic part 5^'' := g^, +5^" . The 
linear part g\. is defined by 



52, fe — -hdyV{^_^ + (^\dgyV(,^_^ + HdyyVl^_^j + \dyyV{^_^ , 



(76) 



qnl t{ 9^^ qnl\ 



whereas, the quadratic part of g^. — {gf ^,g2 k) i^ given by 



qnl ,_ v^ / p. f 

9i,k — 2^ ^i,tCdVij 

(i,j)&J(M.k) 
+ 



Yl ^L ^y^i. 



{i,j)eJ{ALk-i) 



J2 {vL idt)' {devim + {dtf {vim devQ j 

{i,j,l)GX{M,k) 

+ E iyL (dtf (dyvi^m + {dti ivi.m dyvQ 

{i,j,l)eX(M,k-l) 



(i,j)eJ(M,fc) {i,j)eJ{M,k-i) 



+ 



E («{. (5t)' (devl^m + {dtf ivl^m dev, 



{i.,j,l)<^I(M.k) 



^ n 



+ E (<« (^*)' (^.^2,,)(o) + (ao' («2,.)(o) 9.<,) 



(ij',Oei(Af,fc-i) 



Lemma 2.7. [Description of g"^^] Assume that M > 2. The functions g\. and g^" depend only on 
the V; and the dlvj{0, •) where j € |0, fc — 1] and I G |0, [^J J . They are respectively homogeneous 
linear and homogeneous quadratic functions of their arguments. 
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Proof of Lemma 2. 7. It suffices to examine the various terms appearing in the sums involved by 
the definition of g^" . 

- The sums based on the symbol J can be dealt by observing that 

[Af>2, fce{A:-l,fc}, {i,i)(^J{J\l~k)\ => i < fc-1 and j < fc- 1. (77) 

- The sums involving the symbol X are of the form T{M, k) or T{M, k — \). Coming back to 
the definition (41), we can easily infer that 

[M>2, A:e{fc-l,fc}, (i,j,0 eX(M,fc)] =^ i < fc - 1 and j < fc - 1 (78) 

as well as / < [|J . 

D 

Due to Lemma 2.7, the expression g^' can be viewed as a source term in system (75). In order 
to guarantee the fast decaying criterion in r, we can proceed as described below. 

Lemma 2.8. [Solving 15 in the case of a fast decay when t tends to +00/ Select functions 

There is a unique expression Vq — '^{Vq ,Vo ) G {H°°)'' {T x M;M)'^ which can be determined in 
function of G through formulas (80) and (81) such that the Cauchy problem 

{BV^G, V^\V'^,V'//,V'^,V'//), V\^^„=\V,'^X^^X^,K'^^) (79) 
has a global solution V belonging to the space {£^) 
Proof of Lemma 2.8. The strategy is again to argue line after line. 

- First line. Just apply the end of Lemma 2.2. 

- Second line. It suffices to take 

/•+00 
V(-) := - / G^"{s,.) dse {H^)" (80) 



in order to recover after integration that V^^^ G (^J")'^ with 

pT /•+OO 

V^I'{t,.)^V^"{.)+ j G'"[s,.)ds = - j G'/'{s,.)ds. 

- Third line. For all m S N, the operator Ts : if™ — > (Tf™)-"- is continuous. It follows that Tg 
sends the functional space £|" into {8g°)^. Concerning V"^-^ , the argument is again Lemma 2.2 
appfied this time with the source term G^-^ - T,V'^^ - T^V'^// e {£s^)^. 

- Fourth line. For all to G N, the operator Tj : if™ — > (iJ™)// is continuous. Therefore, we 
know that T/ : £g° -^ (fj^)//. With this in mind, it suffices to select 



+00 




^o'^''(-) := - / (G'/Z - TfG^^){s, .) ds e {H°°y . (81) 

D 
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2.1.4 Proof of Proposition 2.1 

The matter is to show by mduction on K E [0, A^ + 1] that the property given at the level of line 
(105) is verified. 

Verification of 'HA/'(0). By convention, we start with v^ = and vl = for k G {—3, —2, —1}. 
Applying Lemmas 2.4, 2.5 and 2.7 with fc = and exploiting the given (linear or quadratic) 
homogeneity properties, we find that /q'- = and gp' = 0. 

Recall that Vq -.— ^^Vq and Vjj := *$Wq. The matter here is to show the existence of functions 
Vo" e CH^)*"* and V"f/' e (£^)*°* such that: 



For M > 3: 



AV,^ = 0, BVj = 0, iVo' + VjM-)^*<^v"oi-). (82) 



For M 



AV.^^O, BV,f=0, iV,^ + V,^){0,-) = '^v",i-). (83) 

By construction, the two first lines of (82) and (83) amount to the same thing as fidgevf-^ = 0. 
We recover here (48). From (68), we can deduce that AVq = AVq. Therefore, the discussion 
concerning (83) is the same as the one related to (82). 

The above initial condition can be decomposed into 

'iVo''^, Vo^'^M •) + \V,^''^, Vj'^M ■)-{!- ^)^^S(•) , (84) 

\vf/^X"^^M ■) + \vi'", Vr"m ■) = Mi-) ■ (85) 

In view of (80) and (81), we must have V^'^'' {Q, ■) = V^'^'' {Q, ■) = Q whatever Vo-''^(0, •) is. It follows 
that we can identify Vq* (0,-) and V^ (0,-) through (85). Now, knowing what is V^'* (0, Oj 
formulas (72) and (74) give access to Vq^^{Q, •) and ^^'^^(0, •). There remains to use the condition 
(84) in order to further extract Vj/ (0, •) and Vj/ (0, •). 

We apply Lemma 2.6 and Lemma 2.8 in the case of the initial data Vq (0, •) and Vj (0, •) which 
have just been computed. Note that, due to the preceding construction, there is no contradiction 
between the expressions Vj'-'-(0, •) and Vq (0, •) thus obtained and the compatibility conditions 
required at the level of (84) and (85). By this way, we can recover functions V^ € (H^)*"* and 
Vq G {£g°y°*. Then, to conclude, it suffices to come back to v^ e H^ and Vq e £g° through the 
action of $^^. 



Assume that the condition HAf{K) is true for some K G |0,7V]. Since the criterion (48) is 
satisfied, the problem can be interpreted as before. The matter is to find two functions V^ i ^ := 
*$w| G (n^Y"^ and Vl^^ := *$w{ G (f^)*"* such that: 
- For A/ > 3 : 

AV^+,^'^ff, Sy/+i=*$5^', (Vl+i +Fi+i)(0, •)=**& <+i(-)- 
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- For M = 2 : 

AV^+i^'^^ff, BVf,^,^'^gl\ (^^+1+14+1 )(0,-) = *'J><+i(-)- 

The induction hypothesis applied with the index K together with Lemmas 2.4, 2.5 and 2.7 say that 
the functions *$/j^' and *$5^' are known source terms with the expected (?^^)*°* and (f^)*°* 
regularities. 

-When M > 3; the initial condition can be decomposed into 

*(i^^Vi-^K4'i)(o,-) + *(^^/+-'i,^#+^)(o,-) = (/-n)<+i(.), (86) 

*(i^k/i' ^^+i)(o, •) + *(^^^^(, v'f^()(o, •) = n<.+i(.) . (87) 

In view of (80) and (81), as clearly indicated in the statement of Lemma 2.8, the expression 
^if+i('-'' ') depends only on *$5^'_,_]^. We can determine V^j_]^(0, •) through (87). Knowing what is 
V^^j^(0, •), formulas (72) and (74) give access to V^+i(0, •). There remains to use the condition 
(86) in order to deduce V^_|_i(0, •) and V^,]^(0, •). Remark that the initial data V^+i(0, •) and 
V^ , j^(0, •) thus obtained inherit the expected H°°(T x M) smoothness. 

Again, we apply Lemmas 2.6 and 2.8 in the case of the initial data V}*/^j^(0, •) and V^_|_]^(0, •) 
which have just been computed. As before, the preceding choices concerning V^_|_]^(0, •) and 
V^ii(0, •) are sufficient to guarantee (87) and (86). We find that V^^^ G {Hff°^ and V[^^ £ 
(£g°y°^. To conclude, it suffices to come back to v^_^^-^^ e H^ and w^_,_]^ G £g° through $^^. 

-When il/ = 2, the same types of arguments prevail. This ends the induction. D 

From the preceding construction, we can also deduce the following information. 

Corollary 2.9. [Nonlinear homogenizationj For all /c G JO, A^ + 1], the expression Hv^ can be 
determined through the following parabolic equation, 

dtiLvi - (a^ + ^n {(d^'h)^)^ dyynv^ - sf , (88) 

where the source term S^' ;= *('S'"fc, •S'j fc) depends on the index j with j < k. This fact may be 
formulated by writing \S^^„ S^,%) - *(/i"i, /2"l)(«g, ■■■,vl-i)- 

2.1.5 Approximated solutions 

Li this subsection, we prove estimate (11). We assume that m > 2 so that if'"(]R^) is an algebra. 

We explicitly compute the action of the operator £° on the approximated solution u° built in 
the previous subsection and estimate the remainder R^ := £"(£,«") in //™-norm. 

To this ends, we justify that the slow profiles u°* and the fast profile v^^ are good approximation 
of operator C"''^ and C^f . 

One aspect of the proof is to compute the difference between £°'''and £"•''* given by some Taylor 
formula. To control it, we provide the following Lemma. 
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Lemma 2.10. Let m > 2 be an integer and S e]0,fi[. Let f G ^^'(T x R), g G -H™' (T 

On the strip [0,T], consider the function h^ [t,-) := f (e^'^t, ■) L u^ g{u, ■)du . Then the family 

{£-(^+i)ft|^p} is bounded m H"'°(T x R), 



xp V 

i.e. 



sup sup 

ee]o,i] te[o,T] 



'^expV^i ') 



_y"(T> 



< +00. 



(89) 



The proof is obvious using some Gagliardo-Niremberg's estimate. 
Proof of Proposition 11. First, we decompose the action of £" on u° into 

£''(£, <(t,.))=/:'^^(e,«r(i,-)) + >C"/(£,«r(i,-),<^(V£'r)), 
where the operators C""^ and C^^ are defined in (32) and (36). 

o Writing w°* as the sum w"* = X]fe=o ^''^fc ^^^ using the cascade of equations (42), we obtain: 

N+l 

k=N 

2(JV+l) + (M-2) 2(Af+l) + (Af-l) 

+ E ^' E <^^K+ E ^' E 

k=N (i,f)eJ{MM+2) k=N (ij)gj-(A/_fc+i) 

/ N+l N+l \ 

k=N k=N 

N+l N+l N+2 

k=N k=N k=N 



^h^V^'3 



K 



(90) 



/ 



First we can factorize by e^ in the above expression. Then, since w^ is in "H^ for fee JO, A^ + 1], 
we estimate (90) in _ff™-norm and get for all integer to > 2: 



sup e sup ||i?^(t, Olljj^jTxR) < +0O- 

ee]0,l] t6[0,T] 

Thus vl is an approximated solution for the operators C"'^ (up to order N). 

o We constructed v^ so that it approximates the operator £"•''* instead of C^^ . To pass from 
C^^ to C"'^'^ , we perform a Taylor formula (with respect to t up to order A^ — 1): 
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where the remainder is defined as: 






t N 
U 



+^ / -^{dtfvf{u,.)dudyvl{t/s^-) 



m 



+ £*^-2 {ti{t/e\-) l^^{d,f-^dovl{u,-)du 



* „,N 



+ vi{t/e\-) j^ —id,r-^dyvliu,-)du 

According to Proposition 1.1, the hypothesis of Lemma 2.10 are satisfied. For all integer m > 2 we 
have: 

sup sup e~^ ||i?f^(i, OlL^fTxR! < +°°- 
ee]o,i] te[o,T] ^ ' 

o Finally, we plug v^ = X]a;=o ^'^''^k ^^^ ''^e — Sfe=o ^'^''^l ™^° ^'^'^* ^^^^ ^^^ ^^^^ cascade of 
equations (43). We obtain: 



rf\e,vi){t/e'')=Rli + Rli. 



The remainder R\' is defined as 



2(JV+l) + (A/-2) 

Rlf ■. = e''dyvf„^,{tle\-)+ E ^' E vl{t/e\-)dovf{t/e^-) 

k=N {i,])eJ(M,k+2) 

2(N+l) + (M-l) 

+ E ^' E <,(t/£^•)5,«/(t/£^•) 

fe=JV (ij)eJ(M+l,fe+2) 

/ Af+1 Af+1 \ 

fe=Af k=N 

N+1 N+1 N+2 

H E e''doevi,it/e\ •) + A E £'5e,«L(t/£^ + A E e''dyyvl,_,{t/e', ■) 

\ k=N k=N k=N I 
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whereas R^^ consists of the terms corresponding to the Taylor formula: 

4N+{M-2) / , 

k=N \{i,j,l)eX{M,k+2) 



{i,j,l)eI(M,k+2) 



k=N \(i,j,l)eI{M,k+l) 



£l(M,k+l) I 



-By construction (see Proposition 1.1), the profiles {ffe}fcg|o,Af+i] lie in £^. Furthermore, we can 
factorize Rl^ by e^. We deduce that for all integer m >2: 

sup sup e"^ ||i?e^(t, OlIff^fTxR) < +°«- 
EG]0,l]te[0:T] *- ' 

-There remains to estimate the term R^^ . A priori this term can be dangerous because it contains 
some polynomials in the variable i/e^. Nevertheless we can use the fast decreasing behaviours of 
the profile v^ . Indeed if / € £^ then for all / e N the function r' / is also quickly decreasing: 
T^ f & £g^ for some < S' < 5. Furthermore, noticing than we can factorize by e^ in R^-^ , we 
obtain for all integer m > 2: 

sup sup e"^ ||i?^^(t,-)||^„fTxBl < +°°- 

Ee]o.i]te[o,T] ^ ' 

D 

2.2 The case of the pressure - Consequence 

Here, we still assume that M > 2. First of all we quickly prove Proposition 1.2. Then we take 
advantage of the control obtained on {q^}e to prove that the approximated solution v^ is a good 
approximation for operator (£^,£^) assuming v is large enough (Proposition 1.3). 



Approximated pressure. Consider the approximated velocity v^ built according to Proposition 
1.1. Since the operator Cq is linear with respect to the pressure variable, we build the profile q^ as 
the solution of the following problem: 

^o(e,gLw£) =dtql+e-^hdyql 

+ e^^-^{vfdeql + e vl-'dyql) +C e^'-^ql{devf + edyvl") = (91) 

with initial data satisfying g|(0, •) — q^{-)- 
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At fixed e, for any positive time T, (91) has an unique solution qf. in H™(T x M). We recover the 
approximated solution q^ (on the strip [0, T]) by summing over all the multi-indices k E JO, iV + 1]: 

N+l 

qt := ^ e\l. 
k=a 

However, since the transport is singular the family of solution {q'^}c is not bounded in ^"(TxM). 
Yet, it is in the anisotropic Sobolev spaces. In other words. Inequality (14) is satisfied. 



Approximated solution for the operator C. A direct consequence of Inequality (14) is that 
{(g",w")}£ is an approximated solution for the operator £, i.e. Proposition 1.3 is satisfied. 

First, we have Co[e,q'^,vf) ~ and: 

' {C^,C2) (£,g^<) = C^ie^v-) + Ce^'^-''-^ * [q'^deq'i^eq'idyq'i) . 

The quantity £°(e,u^) can be estimated thanks to (11). The pressure term e^^^^^^^ *(<7°9e(7",eg" 
dyq"^) can be estimated thanks to the Gagliardo-Niremberg's Inequality. Let a e N^, \a\ < m, then 



Wd'^ (q^deqmL^ = d"do{q^ 



^2 — 9 W'ie IIL°° liyellj^^i^i^i 



To recover an estimate in the anisotropic version of the Sobolev spaces, we use (5) together with 
the following equivalence of norms: 



Lff" 



<£" 



H" 



(92) 



Thus, we obtain 



^-N ||g2iy-A/-2 



{q^dgq^,eq^(Jyq^)\\^„^^e ll'J'e ll//^"^'+J ^ £ 



2i/-M-5/2-(m+l)-JV 



Assuming (15), it completes the proof. ■ 

3 Energy estimates 

In this section we prove Theorem 1.4. The integers tti, ly, M, N and R satisfy property (18). As 
ever mentioned at the level of the introduction, page 12, the main issue is to get a control over the 
singular term £~'^dghv^^ . To desingularize it, we consider the new unknowns: 



rifl 



..IR 



^R 



,.2-R 



It satisfies a hyperbolic-parabolic system (singular in e) : 

dtqf + e-'h dyq^ + s''-' {vl^ dgq^ + £«f 9,§f ) + C e''-'qf' {devl'^ + 8 9,«f) 



QT.R 



e h d,,v, 



V"e 



(0,e ^dghv, 



1R\ 



ev'e'-dyv^) 






(93) 

(94a) 
(94b) 
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together with the mitial data (<;^(0, ■),v^{Q, •)) = 0. The right-hand side of Equation (94) is defined 



as: 



^Ce"+M-'qR^devl^ + d^^^R)^ (95) 



(9,,e2a,)(q-+£V) ■ (96) 

This is clearly a non-linear problem and the variable of pressure qf and the variables of velocity 
v^ are coupled. The dissipation is turned into 



We clearly desingularize the hyperbolic part to a cost on the parabolic part. Q^. It can no 
longer satisfy an estimate such as (28) necessary to keep control over singular terms in (94). One 
important aspect of the proof is to obtain Inequality (108). If it is not satisfied, the mechanism of 
the proof fails. 

Thus, in this section, we look for accurate energy estimates for the unknown {qf^, vf). We easily 
go back to the initial variables according to (93). 



Coupling and nonlinear aspects. In Equation (94b), the variables of pressure and velocity are 
only coupled threw the term: 

"^^ \de.e^dy){qt + e'^qff. (98) 

If V is large enough, the Equation (94b) is somehow independent of the pressure. It seems that we 
can deal with the velocity and then deal with the pressure. This is reinforced by the fact that the 
pressure and the velocity are estimated in different Sobolev Spaces. 

Yet, the term (98) has to be estimated carefully. We have to link the anisotropic Sobolev norms 
with the classical Sobolev norms (92) . It explains the loss of precision on Wm with respect to the 
regularity m (see Equation (18)). 

In practice, we perform energy estimates on Equation (94b). Then we plug the estimates ob- 
tained on the velocity into Equation (94a). We underline here the fact that we need more regularity 
on the velocity to estimate the term S^'^'^ (defined by Equation (95)). Thus, the regularization 
of the velocity (thanks to the dissipation) plays again a crucial role. 

Finally, the problem is nonlinear (for example term (98)). Classically, nonlinear terms are 
estimated thanks to the Gagliardo-Niremberg's estimate which required a L°° control over the 
unknowns. To get this control, we introduce the characterized time T*: 

T; :=minh,^sup^ JvtG [0,T], ||§f (t, .)L..+3(^,«) < 2, ||^f (i, •)L„+3(txk) < 2}! ■ 
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It provides L°°-estimates on the strip [0, T*], for all e e]0, 1] (see (5)): 

Vt e [0, r;], \\V~eqf{t. •)L^-^+.,^(TxE) < 2, \\v^{t, OL^.i.oofTxM) < 2- (99) 

3.0.1 Results and Consequences 

As indicated, we first prove an estimate over the velocity on the interval [0, T*], thanks to Equation 
(94b): 

Proposition 3.1. Let v^ be a solution of (94) on [0,Tg*]. There exist a positive constant Ecrit and 
two positive constants K^ and K^-^ (independent of e) such that 

Vee]o,£,„,], vte[o,r;], ||^f(i,-)||J,„+3(TxK) < ^'""^m (e"^"* - 1) • (loo) 

Furthermore, one can prove a regularization property. Select a multi-index a G N^ of length \a\ 
smaller than m + A; then 

V£e]0,e,„t], Vte[0,T*], / e-(i-*°i-«)9"^f(s,-) ' ds < K^^s^'^'-t , (101) 

where Si,j denotes the Kronecker symbol (of two integers): Sij — ii i ^ j and di^i — 1. 

By plugging (lOO)-(lOl) in (94a) we obtain: 

Proposition 3.2. Let q^ be a solution of Equation (94) on [Q,T*]. There exist a positive constant 
Ecrit o.iT'd a positive constant K^ (independent of e) such that 



VeG]0,£c„*], Vte[0,T;], \\q^it,-)\\l^^, <E^^-K, 

(1.=) ^ ' 



t. 



Those two lemmas allow to prove an accurate estimate of the solution on a strip independent 

of E. 

Corollary 3.3. Consider integers m, u, M, N and R satisfying the condition (18). Then, there 
exist two positive constants Ec and Tc (independent of e), and a positive constant c^rr > such 
that, 

VigpjTc], V£e]0,£c], Pe (i, Oll^m+S^TxK) - '^err, || W^ (i, •) ||^„+3|.^^jj^ < Cerr • 

Proof of Corollary 3.3. We argue by contradiction: 

V(e,r)e]o,i]x[o,i], 3£e]o,?], t; < r. 

We recall that there exist Ed and C^ positive constants such that, for all e e]0,erf] and for all 
time te [0,r;]: 



\n^(i .\\\ ^ 2wm r-i 4. llTr^rf ."ill <^ 2to„^ . 



We choose for instance T = min( ^p , 2) and e = e^ < 1. In particular, T < 1. From assumption, 
there exists Eq ^]0,Ed] such that T*^^ < T. Furthermore, since Wm > from condition (18) we get: 

vtG [o,r;], ||C(^'-)L™+3^(Tx«) < 2 <2' ll^fo(^'-)L.+3(Tx«) < 2 <2- 
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te[0,T;j^||gf„(t,.)L-+3 ,,^) and i e [0,T;J^ ||^fjt, .)L„ 



3(TxR) 



Now consider the applications 

They are continuous, qf' (respectively v^) can be extended in time as long as the quantity g^ (t, •) „„i+3 

( respectively wj?(t, •) Lfm+s) remains bounded. 

It follows that we can find T £]T*^,T^^\ such that for all time t e [Q,T], ||g^(t, ■)||h"+3 < ^ 
(respectively llu^ (f , .) 11 „„+3 < 2). This is in contradiction with the definition of T* . D 

Assuming Propositions 3.1 and 3.2 are satisfied, Corollary 3.3 holds. We go back to the initial 
unknowns thanks to the Equation (93) and obtain Theorem 1.4. 

Thus in Subsection 3.1, we prove Proposition 3.1. We take advantage of Inequalities (100) 
and (101) to obtain Proposition 3.2. 

3.0.2 Notations 

Here we introduce some notations required for the proof of Propositions 3.1 and 3.2. 

Let (q^, u^) an approximated solution of order TV constructed on the interval [0, 1] according to 



Proposition 1.1 and Proposition 1.2. The family {(g^jV")}^ lies in 'H^7[ \ x "H™ ' ^^id satisfies 
Inequalities (11) and (14). We denote by Ca, Cc^ and Cc positive constants such that: 

sup sup ||w°(t,-)||ff™+6 <Ca, sup sup ||g,"(t,-)||a™+6 <Ca, (102) 

ee]04] te[o,i] e6]0,l] te[o,i] <^'=^ 

and 

sup sup ||e-^/:'^(£,0||^„+e <C£. , (103) 

ee]o,i] te[o,i] 

sup sup We-'^ £o{£,qe,v^)\\H^+e„^„.<Cc. (104) 

ee]o,i] te[o,i] <^'=' 

In what follows we adopt the conventions: 

VmeN*, ||/||^_=0, 

to simplify all the statement. 

3.1 Energy estimates for the velocity 

In this section, we prove Proposition 3.1 by induction on the size of m settings: 

V{m) : " Proposition 3.1 holds up to the integer m". (105) 

To go from ?7^ to m + 1, we prove an energy inequality for the velocity performing an energy method 

o 

on Equation (94b) in the homogeneous Sobolev space iJ™ (defined page 3). 

36 



Lemma 3.4. There exist ci and Ed two positive constants (which only depend on m) such that for 
any J £ |0, r7i + 3], there exist four positive constants Cp, Cj, Cj and Cj such that for all e g]0, Ed], 
and for all time t € [0, T*], 



+ "-"J Ire l^^J Vll^j-i(TxR) + '-'./ Z^ II ^eV^^ J 



L2(TxI 



Subsections 3.1.1-. . .-3.1.3 are dedicated to the proof of Lemma 3.4. Then in Subsection 3.1.4, 
we finally prove Proposition 3.1. 

First, consider a multi-index a e N^ and differentiate a times Equation (946). Then we multiply 
by d'^v^ and integrate (with respect to the space variables 9 and y), 

]-dt ||a"vf 11^, - (9"(Qewf ),5"wf ) = {d°'Sf'^,d"vf) - (9" {Avf),d°'vf) 

- (d" (Swf), 9"vf ) - (9" (Cwf ), a"wf ) - (a" {Hvf),d°'vf) , (106) 

where the operators T-L, A, B and C are defined as follows: 

Uvf := e-'hdyv^, Av^ := e''~^ {vl^de + i^^9,) '(«]", £t^f), 

Cv^ :=* {0,e~^dehvl'') , Bvf:= e''-^ [vl'^devf + evl'^dyvf) . 

The strategy of the estimates is the following. We first prove the coercive estimate over Q^ (see 
Lemma 3.5 in Subsection 3.1.1). Then we compute all the singular contribution (with respect to 
the regularity and to e) in Lemmas 3.6-. . .-3.8 that we absorb thanks to Q^ (see Subsection 3.1.2). 
Finally, we simply bound the remaining terms thanks to Lemma 3.9 in Subsection 3.1.3. 

3.1.1 Step one : Coercive estimate for the dissipation 

We start by estimating the term involving the dissipation Q^. 

Lemma 3.5. We recall that /i > 0. Select constant A and fj, which satisfy 

A<4^. (107) 

There exists cg > 0, there exists Ed g]0, 1], such that for any function f G H^{T x M,M^); 

V£e]0,£rf], -(Qs/,/)>Co*e(V,/). (108) 

where $j is defined in Equation (28). 

Proof of Lemma 3.5. Let / G i/^(T x M,M^). We decompose / in Fourier series (in 0) into: 



m2/):=E/fc(2^) 



^Me^ 



k£l 
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then we use a Fourier transform in the variable y. By the Parseval equaHty, up to a constant (that 
we forget here), 

- (Qe/,/) = E / M)Qe{KOM)d^. M) := / e^y-^h{v)dy, 



where Q^ is defined as: 






We interpret {Qef,f) as a quadratic form in the variables /^ and /|. In this way, to prove Inequality 
(108) we show that there exist Ed and cq two positive constants such that for all g G iJ^(]R,IR^): 

VfceZ, VeeM, 5(OQe(^,05(C)>co(e-'A;2+e2) (52(^)+52(^)) , (IO9) 

At fixed (fc,^) e Z X M, Q^{k,£^) is a diagonalizable matrix since real and symmetric. We compute 
its eigenvalues: 

^^l{k,0 :- ^^^^{{e-'kr+e) - y V'((^"'^)' + ^') + (^(e^'fc)^^ + ^-^ [e-'k^^ • 
In the end of the proof, we show that for any (fc, ^) G Z x M we have 

^Jil{k,0 > Co {e-'k^ + e) {91(0+ dliO) ■ 

Since /i^ > /^?, we clearly obtain Inequality (109). 
o We define the function /i^ : R'^ — > R by: 



^ie{x,y) := ^^ "^ {x^+y^) - y v^(.t4 + y4) + (ea;2;)2 +e-2 (xy)^ 
Then there exist e^ and cq (independent of e) such that 

V(.T,2;)eM^ f^e{x,y)>co{x^ + y^). (110) 

-First of all the function /i^ is homogeneous of order 2 in the sense that it satisfies for all a e M: 
V(a;,y)eMxM, fi^{ax,ay) ~ a^iJ^{x,y) . 

Thus, we prove (110) on the restricted set {x, y) G §^ the sphere of center and radius 1. 
-We expand ^.^ as e goes to 0+: 

Me(2:,y) = (a;2 +y2) _ _\x\\y\ + 0(e), 

where 0(e) is uniform in {x,y) £ S^. Let us recall that we have: 

V(a;,y)eM^ (a;^ + y^) - c|a:||y| > <^=^ c < 2. 
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That is the case if and only if A < 4/i, i.e. assumption (107) is satisfied. We get an uniform bound 
of lie (in £)• Finally, there exist Ed and cq two positive constants such that Inequality (110) is 
satisfied. 

o Plugging X = e^^k and y — (, in Inequality (110), we obtain for all e £]0,eci]: 

V (fc, C) e Z X M , fi^k, > CO {{e-^kf + e)- 

We deduce that 

fkiOQAk,OfkiO>coiie-'k)'+e)\\fkm\'- 
Finally, summing over k d 1, and integrating with respect to ^, we have: 

-{QefJ)>coY.f{{e-'kr + e)\\M)rdC = co'i>e{^,f). 

7, ^'2' " M. 



fce 



n 



Remark 1. The singular change of unknowns desingularizes the hyperbolic part of system (94) to 
a cost on the parabolic part Ve- Qe is dissipative enough assuming the A-part of the dissipation is 
strong enough with respect to the V div-part (X < 4/i/ 

Remark 2. The Inequality (108) can be proved for a more general family of dissipation. We can 
replace eX in the dissipation Q^ by some coefficient X^ going to when e approaches 0. Then 
replacing assumption (107) by: 

limsup — < 4/i. 

e-)-0+ £ 

Inequality (108) still holds. 

From Equation (106) and the preceding lemma we deduce that for all e €]0, Sd] and for all time 

ie[o,r;]: 

- (9" {Bv^), a"vf ) - (9" (Cvf ), a"wf ) - (9" {Hvf),d°'vf) . (Ill) 

3.1.2 Step two : estimates over the singular terms 

We now try to absorb the singular terms in (lOG), thanks to $j. Singular terms are of two types: 
-the operator only singular with respect to the regularity, S^'^ , 
-the operators at least singular with respect of e, "H and C. 
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Contribution of S^''^ . We start by estimating the term S"^'^. It contains nonlinear terms 
together with the couphng terms (98). What can be underlined is that this term (98) is nonlinear 
(only) singular with respect to the number of derivatives acting on the pressure. This is a problem 
since the dissipation only regularizes the velocity. When possible, we pass those extra-derivatives 
onto the velocity (thanks to an integration by parts). One can prove: 

Lemma 3.6. [Control over S^'^ ] Select a multi-index a G N^ with length smaller than m+ 3 and 
a positive constant Cs- There exist Cg, Cp two positive constants such that for any £ G ]0, 1], for 
any time t £ [0,Tg*], 

|(a"5f -, a%f )| (t, .) < Cs \\v^\t^ OllJ^M^.f.xM) + ^s \mt^ OIILktxm) + ^^ ''""^ ■ 

Proof of Lemma 3.6. Consider a G N^ satisfying \a\ < m + 3. We decompose the source into 



£ 



R+M-2 ^ga (~lfl g^~R ^ ^ ^fl g^-R^ ^ ga~R^ 



- C2-l£2''-*^-^-2 (t^g^^^g^^ga f^^a ^ ^R ~R^^ ^ ga~R 

We estimate the three above contributions. 

o The first term: contribution of the approximated solution. u° is an approximated solution 
for the operator £". According to Proposition 1.1, it satisfies Inequality (11). That is to say. 
Inequality (103) holds and: 



IgAf-fl 



(a- (£-~£'^(£,<)) ,9"^f)| < I (£2(^-«) ||9" (£-^/:'^(£,<))||^, + \\d"v^\Q , 



< 

o The second term is the contribution corresponding to the non linear part "u • Vw" (in Navier- 
Stokes Equations). Select a positive constant ci. Further, we choose it so that the contribution 
of Wv^ in iJ™-norm is small with respect to $g. 



£ 



R+M-2 I ^^a (^Iflg^^fl + £ ^fl g^yR)^ ga~R^ | 



r2 J 



< ci£2(«+^^-2) (||a" {vl^dov^)\\l, + ||9" {v',^dyV^)\Q + ^ p'^vml, . (113) 

We deal with the nonlinear term d"{vl^dev^) applying the Gagliardo-Niremberg's estimate to- 
gether with (99). There exists a positive constant Cg which only depends on m such that for all 
£ g]0, 1] and for all time t G [0, T;]: 



, 2 



\0 [v^ dev^ )\\^^ <Cg \^\\v^ IIl~ lr«"£ ll^i^i + 11^6 ll^i„i lr«"e I 

< 4c2 fulfil o + ii^f II o Y < 8c2 fulfill, ,,, + ll^f Hi ,) 
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The same holds for w|^ dyV^. Finally, we get that for all e e]0, 1] and for all time t E [0, T*] 



< 16ci C| \ml,,,, + (16ci C'y + ^) ll^f |||,„, . (114) 

o The last term. This term has to be studied with care. As mentioned, two obstacles have to 
be overcome: 

-It is singular with respect to the number of derivatives acting on q. However the term can 
formally be written under the form: 



1 



(V(g^),«) 



A derivative can be passed onto the velocity with an integration by parts. 

-The pressure is only estimated in anisotropic Sobolev spaces HJ^ ^^ . We pass to the anisotropic 
Sobolev spaces HV^ s thanks to Inequality (92) at a cost on the precision Wm (see Equa- 
tion (18)). 

First we integrate by parts: 

We select a positive constant C2, then 

/-i2-lp2i/-Af-J?,-2 It 



< 



(^{de,edy)d'^{q^,+e''qf)\d"vf) 

n2 p2{2}j- R- M -2) 



8C2 



a"(g,"+£V) „ + tII^"*(^«'^^>. 



L2 



e llL2 



(115) 



We apply the Gagliardo-Niremberg inequality together with the equivalence of norms (92) , 



9"(g,°+£V) 



..a I ^R~R- 1 



L2 



< on \\n°- -I- e"-V 



ke+e""??! 



H\' 



< 2Co£-l"l llg" + e^gfll ,^ ||g° + e^of 



e \\L° 



9e ||ol° 
H, 



With regards to the construction of the characterized time T* , (99) and (102) holds. It results in: 

Vse]0,l], VtG[0,r;], ||Vi9f(i,-)|L„.(TxE) <2, ||Vig,"(t,-)|L.(TxK) < Ca- 

We obtain 



d'' {q^ + e'^q^y < 2C<,£-i/2-l"l(Ca + 2)2 



(116) 
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Finally plugging Inequality (116) in (115), we deduce: 



< 



(J2 (j2 g2(2i/-_R-M-5/2-|a|) 



2C2 



(Ca+2)4+|||.,^ll^^ 



ff|a| + l 



(117) 



o To finish, we put estimates (112), (114) and (117) together. Let ci and C2 be two positive 
constants for all e g]0, 1] and for all time t e [0, T*]: 



|(a"5f ^,9%f )| < (I + 16ci cl) ||tf |||m.^ + (^ + 



16c, C| + ^ 



„-R| 



l_H-| = l 



+ (Cc^ e^CW-fl) ^ ^^(C, + 2)4£2(2V-fl-A/-5/2-|a|) 



<(| + 16ciC^2 



lffl°l+l 



2+^^^^^^^+2^ 



"e WhM 



+ (^C.. + ^(a + 2)^je^ 

Select C5 a positive constant. We choose ci and C2 two positive constants such that C5 := ^ 
16 Cl C|. Then, it requires C^ = ^ + 16ci C^ + ^i^ and Cp = C^a + ^{Ca + 2)^. 
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Contribution of operators "H and C. Presently, we study the singular operators Ti. and C. We 
decompose their action into: 

(a" (Cwf) , c'^iyf ) = (c (a"wf ),^f ) + ([a",c]wf,c'"iyf ) . 

by making some commutators appear. Since the commutator of two operators respectively of order 
m and n is of order m-\-n—l, the contribution of the commutators are less singular than expected: 

Lemma 3.7. Select a multi-index a G N'^ with length satisfying 1 < |a| < ttt, + 3. Select V G 
{HjC}. There exist Cy, Cy and Cy three positive constants such that for all e e]0, 1], for all time 

te[o,r;], 



e'i^'S''i.o)d°'vf{t,-) 



Cl 



2 

L2(TxI[ 



The proof is rather easy. We do not write it. Then, since {% d^v^ , vf") = 0, there remains to 
deal with {C{d°'v^),d"v^). 

Lemma 3.8. [Absorption of C] Select a positive constant c\ and a multi-index a G N^ of length 
less than ttj + 3. Then for all e G]0, 1] for all time t G [0, T*\: 



(ca"^f(i,.),a"iyf(t,.)) < ^%^l|9"i^f(<r)lli.(TxM) + | 11/^11 



^oo ||£ ObO V^ 11^2 



42 



Proof of Lemma 3.8. We prove the result for a = (0,0). Replacing vf' by d°'v^ in the above 
estimates proves the result. We select ci a positive constant. The idea is to integrate by parts with 
respect to the variable 9 to make the weighted derivative of the velocity e^^dgv^ appear. 



\{Cv^,v^)\^ 



I he-^devl^vl^dOdy + j hv]^ e'^dg'i^^^de dy 



< 



\h\\ 



|li^||^2 +Ci ||£ OffW 



2 
e llL2 



D 

Absorption of singular terms. Select ci and Cg two positive constants to be chosen (small) 
later. We get a bound for the right-hand-side of (111) applying Lemmas 3.6-. . .-3.8. We obtain 
that for all e e]0, £d] and for all time t e [0, T*]: 



i9,||a"^f(i,.)ir^.+ 



co$e(V,9"^f)-f ||/i|| 



2 ^ \\~Ru \\\'^ 

e llL2-Cs||V(^,-)||^,+i„ 



-iafl9"w^ll -r<^\\r,R 



(Txl 



< 






2w„ 



+ (Q' + c^) ll^f (t, ■)||h..-(txe) + (^c + cl,) £-(i-^-«)a"^f (i, •) 



L2(T> 



-h I (a" (y^vf ) , 9"wf ) 1 + 1(9" (Svf ) , a"wf) I . (118) 



At this stage, the quadratic form $^ can not absorb the norm ^^^(i, •) ° • However sum- 

*= ' ^ '^ II e ^ ' ^lli/J + i(TxR) 

ming over all a G N'^ of length J the quadratic form X]|q|=j *^£(^; d°'vf) does. 

Choosing ci and C^ such that the quantity cq ;= cq — ^ II^IIl°° ~ ("^ + 4)C5 is positive, we 
have: 



^ (co$, (V,9"^f) - I ||/i|l^^ \\e-^ded'^vff^,-Cs \\v^{t,-)fl,,,. 



W\=J 



(TxR) 



>co ^ <i>,(v,9"^f)-|l|/i||^^ ^ ||£-ia,a"«f||^,-(j + i)C5||«f(t,-)| 



|a| = J 
Cl 



\a\ = .J 



ff'^+i(T> 



> (co-y||/i|Lo.-(m + 4)C5 ^ $,(V, 9"^f ) = ?o ^ 'i>e(V, 9"^f ). 



|a| = J 



|a|=J 
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Hence, we obtain for all e GjOje^] and for all time t E [0,T*]: 

|a|=,7 



2 

L2(TxI 



+ J2 \{9°'{^v^),d"vf)\+ ^ |(9"(Svf),9"wf)|. (119) 

Remark 3. The term X]|q|=j k~ ^~ ^^'"^"^if (^i ') r2fT lu'i ^^^^^^ ^^ ^^ singular with respect to 
e. With regards to the Inequality (101), we will interpret it as a non-degenerate source term in the 
proof of Proposition (3.1). 

3.1.3 Last step: estimates over non-singular terms 

The two last terms (contribution of operator A and operator B) can be easily estimated. They 
satisfy: 

Lemma 3.9. Select a multi-index a G N^ of length smaller than m + 'S. Select V G {yl, S}. There 
exist Cy and Cy two positive constants such that for all e g]0, 1] and for all time t G [OjTg*]; 

1(9" (V^f) ,9"^f)| (t) < Cl ||^f(V)|||M(,,„) +C$ ll^f (V)lliM-.(TxE) ■ (120) 

The proof is very classical. We do not write it. Applying Lemma 3.9, we obtain from Inequal- 
ity (119) that for all e G]0,£d] and for all time t G [Q,T*]: 

^^*ll^e^(V)|lL(Tx«)+^o E *4V,5"i^f) <C}||iyf(V)||L(,,,)+C,(J+l)e-" 

|a|=J 



L2(TxI 



With C] := (J+1) (\^^+Cl + C\ + Cl + Cl+Cl,), C^j := [J + l){C\ + Cl + Cl + Cl) and 
That achieves the proof of Lemma 3.4. ^ 



3.1.4 Proof of Proposition 3.1 

This subsection is dedicated to the proof of Proposition 3.1. We prove by induction that prop- 
erty 7^(J), defined at the level of (105), is satisfied for J G |0,to -f 3]. We proceed in two steps. 
First we obtain an estimate for the velocity (Inequality (100)). Then, we take advantage of the 
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quadratic form J2\a\=j^£ (V,9"w^) to obtain the regularization Inequality (101). 

Since the proof is exactly the same for the case J = and the increment of the induction, we 
do not write the details to prove that property V{0) is true. 



We assume that 7^(7) is true for J E [0, r7i + 2]. We apply Lemma .3.4 in the case J+ 1 (> 1). 
There exist Ed and ci two positive constants, there exist C}^]^, Cj_^^^, Cj^^ and Cp four positive 
constants such that for all e E]0,ed] and for all time t G [0,Tj,*]: 



ldtm^r)\\U.„^,,+c, E $e(V,a"^f)(t) 



|a| = ./+l 



^ '^J+1 |re'(^' ')||/f.7(TxR) +^J+1 Fe (*' OH^j+i 






L2(TxI 



(TxR) 

+ (J + 2)Cp£2-™. (122) 



o First, we look for the iJ''+^-estimate. We neglect J2\a\=j+i 'i'elV, (?"f^) since it is positive 
as a sum of positive quadratic forms. Hence, for all e E]0,eci] and for all time t G [0,Tj?]: 



where Sj-^-i is defined as: 



Sj+i{t) :— Cj_|^i \\v^ [t, •)||i/j(TxR) 

|a|=J+l 



L2(TxR) 



+ ( J + 2) Cp e^ 



From the assumption V{J), Inequalities (100) and (101) hold. The functions t i — > ||?;^(i, •)||^ 

are in L^{[0,T*]). In addition there exists M^ a 



andt^^ J2 £"^^"''°i'"^a°wf(t,-) 

positive constant such that for all e G]0,ec;] and for all time t G [0,T^*]: 



L2(TxI 



1 ^2u)„ 



v^itr) rr,<M^e 



and E ||e"*^^'^"i'''^5"vf(t,-) 



L2(T> 



< ML e 



1 ^2w„ 



Thus, the function t i — > Sj+i{t) is in L^{[0,T*]) and satisfies for all e G]0,£d] and for all time 

tG[o,r;]: 

Sj+iit) < {2Ml + (J + 2)Cp)e2-™. 
We can apply the Gronwall's lemma, for all e G]0,ec;] and for all time t G [G,T^*], 

ll^?(*'-)lH...(TxR) ^ /e^^5..(-«)5,+,(.)ds < '^'"'+j/_+^^^^ e^-- (e^^.'-* - l) . (123) 
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Setting K]^-^ := (2M^ + ( J + 2)Cp)/(2C3_^i) and K],^^ := 2C3+i, the first inequality of V{J + 1) 
is proved. 

o Then, to obtain the estimation i^([0, T^*], i7'^+^), we go back to Equation (122). We integrate 
it with respect to the time t, 



•^ I „, I Til 



|a| = J+l 



<^J+1 



i/ (ll^f(5,-)||^.+,+^,7+i(s))rfs. (124) 



Since ||wf^(i, .)|| o j+^ is positive, we neglect it in the left-hand side of the Inequality (124). According 



^H 



to (123), the functions 1 1 — > ||w^(i, •)|| ° j+i and S'j+i are in L^([0, T*]). Furthermore, there exists 
K'j^^ such that for all e e]0,ed] and for all time t € [Q,T*]: 

\\vf{t, .)|||,,, + Sit) < K]^,e^-- t < Kj^.e'--. 
So to say, there exists Kj^-j^{> 0) such that for all e g]0, Ed] and for all time t € [0, T*]: 



Jo ,. , , , 



'° \a\ = J+l 

We end the proof choosing Kl^ :— max ii'* (i £ {1, 2}). D 

jel0,m+4] -^ 

Remark 4. /n the above discussion, we only use the assumption M > 2 whereas we assume M > 7/2. 
It becom,es crucial when estimating the pressure. 

3.2 Control over the pressure 

In Subsection 3.2 we prove Proposition 3.2. The result is again proved by induction on the size m 
setting Q{m,): 

Q{m) := " The Proposition 3.2 is satisfied up to the integer m". (125) 

To go from rn to m + 1, the proof is once more based on an energy method for Equation (94a) in 
the anisotropic Sobolev spaces H?l ^\ (defined page 4): 

Lemma 3.10. There exits a positive constant e^ > such that for any J G JO, tti + 3], there exist 
three positive constants C], C^j and C^ such that for all e ejO,^^;] and for all time t G [0,7^^*]; 

+ C5(||iyf(t,.)||U^(,,„)+£^(^-^)). (126) 
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The lack of dissipation has two consequences. 

-There is no absorption phenomena. It makes us consider the anisotropic Sobolev spaces instead 
of the classical Sobolev spaces. An unexpected effect of considering the anisotropic Sobolev 
norms is that the family {deq^}e becomes singular with respect to e in L^. However {edeqf}e 
is. We introduce a power of e when necessary thanks to the integer M. It explains why M is 
assumed at least larger than 3. 

-There is no regularization phenomena over the pressure. A difficulty appears when we want to 
estimate terms such as Vp. Each time it appears, we integrate by parts (if possible) to pass 
the derivative on the velocity. It explains the appearance of the H™+^ norm of the velocity. 
The regularization of the velocity plays again a crucial role in this process. 

The energy method in HV^ . consists in differentiating Equation (94b) by e°'^d°'. Then we multiply 
it by e"i9"5? and integrate (with respect to the space variables {d,y)), 



1 9 



-^dt \\e°"d°'q^\\l, + (e"ia"CHgf ),£"i9"gf ) 



= (e"i9"S'^"^-«'^,£"i5"qf ) - (e"i5"(JTJf ),e"i9"gf ) , (127) 
where operators %, B, D and J- are defined as 

Uqf := s-^hdyqf, Fvf := e^'-^ [v]^ deqt + s ??^ 9,9,") + C e'^'-\'i {devl^ + e 9,S?«) , 



Bq^ := e*^-2 (^la Q^~R^^ ^2a Q^~R^ ^~R .^ ^ ^M-2~« ^g^^la ^ ^ g^^ 



2a\ 



In the sequel we prove several lemmas where we estimate each contributions. We sort them 
contingent on how much they are singular with respect to e in the anisotropic Sobolev spaces. In 
Subsection 3.2.1, we study the contribution of operators H, B and V. Then in Subsection 3.2.2, due 
to technical computations, we have to assume that Af > 7/2 and deal with the terms J^ and 5"°'^'^. 
Finally in Subsection 3.2.3, once Lemma 3.10 is proved, we present the induction to prove Propo- 
sition 3.2. 

3.2.1 Non-singular terms: M > 3 

First we deal with the contributions of operators H, B and V. Since we introduce the anisotropic 
Sobolev spaces to get ride of the singular transport H; it is no longer singular: 

Lemma 3.11. Assume that M > 2. Select a multi-index a G N^ of length smaller than m + 3. 
Then, there exist two positive constants C^ and C|^ such that for all £ g]0, 1] and for all time 

ie[o,r*].- 

When \a\ = 0, the contribution of "H is even vanishing. Since it is very classical, we do not write 
the proof. 

Presently we compute the contribution of operators B and V. We have to estimate the term dgq^ 
singular in e in the anisotropic spaces. It thus requires: 
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Lemma 3.12. Assume M > 3. Select a multi-index a G N^ with length smaller than m + 'S. Select 
V G {B^D}. There exist two positive constants Cy and Cy such that for all e g]0, 1] and for all 
time t G [0,T*], 

|(£-5" (V^ ,e-9"gf )| (t) < C^ \\q^{t, .)|||,„, + C^ ||§f (t, •)||i,=i-(,,„) • 

Of course the contribution of D only requires Af > 2 since dyq^ is bounded in the anisotropic 
Sobolev spaces. 

Proof of Lemma 3.12. We prove the result for operator B. We consider a multi-index a G N^ such 
that \a\ > 1. We have: 

The two terms are estimated performing the same proof. Thus we only write the proof for the term 
gAf-2 ^gQi^a ^^ai g^—fl-j gQi^a-fl^y Of coursc, the estimates for the second term only requires 
M > 2 since it appears the derivative dy instead of dg. 



We apply the Leibniz formula. Then to diminish the number of derivatives acting on the pressure 
(for the extremal term), we consider an integration by parts. There exists a family {Ca^/s} of positive 
constants such that, 

gAf-2 l^^ai^a (vfdeq^),e°''d°'q^)\ 






_^(Af-2)/2/ devf{e'''d''q^fdedy. (128) 

JtxS. 

o First, we get a bound for the last term in the right-hand-side of Equation (128). For all 
e g]0, 1] and for all time t G [0, T*]: 



_g(A/-2)/2 / ae«f (i,-)(e"'5"qf(t,-))'d0dy 

JlxR 



<ll<(i,-)IUi,oo/2||g,(i,. 



ff ' 



<Ca/2\\q,{t, 



•h' 



(129) 



o Then, to control the other terms in Equation (128), we make appear the e-derivative edg 
paying a loss of precision on M . For all e g]0, 1] and for all time t G [0, T^*], 



.A/- 2 



-2ai 



a"-^ff (t, •) d^deqfit, •) a"gf (i, ■)dedy 



Txl 



-A/-3 



.2qi + 1 QQ-/3 al 



d'^-^vfit, ■) d^dgqfit, ■) d^qf{t, ■) dO dy 



< 



HHtr)\ 



W"^+'>"' 



+ia^aegf(t,-)|l^ + |k"^9"gf(i,-) 



< 



Ca 



f||5f(i,-)||^i = i-i+2||§f(i,-)||^i»i ) 



Choosing A := J2B<a C^a,pCa/2 > and B := {J2B<a ^a,i3 + ^/'2)Ca > we obtain that for all 
£ e]0, 1] and for all time t G [0, T;], 



|2 . „ii~p,, oi2 



That ends the proof. D 

3.2.2 Singular terms: M > 7/2 

There remains to estimate the contributions of T and 5^'^'^. They are put aside since we 
need M > 7/2 to desingularize it. It is actually a technical assumption and it should be re- 
laxed to M > 3 by proving more regularity on the approximated solution (g^,ti°). 



Estimate for operator T. We start by estimating the contribution of the operator T. The 
main difference with the three previous operators ("H, B and V) is that the operator J^ acts on the 
velocity v^. We have to be prudent with terms of the form 

First it is more singular with respect to the number of derivatives acting on the velocity. Once 
more, when estimating the product, we have to use a control over the pressure q" in L°°-norm. 
We may lose some regularity to go back to the anisotropic Sobolev norms. This term can be more 
singular in L^-norm. Thus it requires: 

Lemma 3.13. Assume M > 7/2. Select a multi-index a G N^ with length smaller than m + 3. 
There exist two positive constants C^ and Sj- such that for all e g]0, 1] and for all time t G [0,T^*].' 

|(e-9" (J^f ) ,£-5"gf )| it) < S^ \\v^{t, •)||L™..(TxM) + G> l|5f(^' OII^M ^(,,«) . 

Proof of Lemma 3.13. Select a multi-index a G N^ satisfying \a\ < m + 3. We decompose J^ as 
follows: 

+ (C£^-2e«ia« {q^devl") ,e"^a"§f ) + {Ce^^~^e"'d" {q^^edyvl^) ,e"'d°'q^) . (131) 

The first and second term (respectively the third and fourth term) can be computed following the 
same line. So we only provide the estimates for the first term (respectively the third term). 

o The first term. We start by estimating the contribution of (^e'^^~^e"^d°' {vl^ deq^) , e°'^d°'qf): 



^^M-2^arQa (^^R Qgq'^) , 8°"' d^'qf) < ]- (\\e'''^''^ e"' 8°" ivl^' deqt)f^2 + ||e"^a"gf | 



49 



We use the Leibniz formula. Then apply the Minkowski inequality and loss a power of e to get a 
control of the approximated solution in the anisotropic Sobolev spaces. There exists a family of 
positive constant {Cq^^} such that 



\^M-2^aic^a f-lR. fi n\\\ ^^M-2 \^ ^ llc-ai + l-l^^^/^ ^o|| W f)Oi-Pz'^R\ 

\e SO [v^ Ceg^Jll^aSe ^ '^a,i3\\£ o'^Oeqi;\\j^^\\o ^v^ | 



L2 ' 



0</9<a 



< 



-J\/-3 



X] C'a:/9 lk£llw(7+''-°= ll"f llff- 



0<;3<a 



Then, we use the imbeddings \\q'^{t^ •)||^m+4.oo < e 2 ||q"(t, •)||^m+6, 

(1,e) (l,e) 



^M-2^ai Qa /~1_R Q ^a^ II ^ ^M-7/2 \^ r' WrflW IU,-R|I 



<a 5] C„,;3||t'f(i,-)L„ 



0<^<Q 



Finally for all £ e]0, 1] and for all time i e [0,Tg*], 



\ \0<,9<a / 



I £ \ ' ''IIh'" + 3 ^ II "e V ' ''I 



L2 



(132) 



o Estimate oj the third term. We now move to the estimate of (C e'^^ ^e°'^d°' (9° dgvl^) , e°'^d°'q^). 
First, apply the Cauchy-Schwarz inequality together with the Young inequality: 



|(C£*^-2£"i9" {qtdevl^) ,£"^a"gf)| < ^ (\\C e''^~^e°''d°' {q^, dgvl"-)f^^ + We'^'d^q 



■R\ 
e I|l2 

To get rid with the contribution of the nonlinear term ||C£*^^^£"i9" (g" Qgw^^) [Lj we start with 
a Leibniz formula, then: 



O<0<a 

< C£*^-2 J2 C^^0\\e°''d'^q^\\^^\\d"-^dgvm^, 

O<0<a 
^ ^ A.f-2 Y^ ^ II an II Rll 

< Le 2_^ Ca,i3 \\q^ 



0</3<Q 






-Af-5/2 



E '^".'9 Il9elli?(7+^ lh?llff'"+-« ' 



< c, 

0<,9<a 
0<P<a 
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Finally for all e e]0, 1] and for all time t E [0, T*] we have, 



2 



< Ic'Cl Yl C^A ll-fllU^ + J||5f|||M . (133) 

\0</3<a y ^ ''■^' 

It proves Lemma 3.13. D 



Estimate of the source term S^'^'^ . The term S^'^'^ contains all difficulties ever met: 

-Singular terms (with respect to the number of derivatives) such as Vgf^ are dealt by integrations 
by parts to a cost on the regularity of the velocity. 

-Singular terms in e in anisotropic spaces such as d$q^ are dealt assume the integer M is large 
enough. 

Lemma 3.14. Assume M > 7/2. Let a multi-index a G N^ with length smaller than m + 3. There 
exist three positive constants Cg, Cg and Ss such that for all e g]0, 1] and for all time t £ [0,T*]: 

+ C's (l + he (*' •)||h-+4(TxR)J ll^e (*' ') ||^M-1(txR) 

Proof of Lemma 3.14- We decompose S^'^''^ into: 



+ (C£^+*^-2^"i9" {q^dgvl") ,e"i9"gf) + (C e^+^''~'^e"'d°' {q^ dy^,^) ,e°''d°'q^) . 

In what follows we study each of these contributions. Since contributions of the second and third 
term (respectively fourth and fifth term) are estimated following the same steps, we only write the 
estimates for the second term (respectively the fourth term). 

o Contribution of the First term. Since q^ is built as an approximated solution it satisfies 
Inequality (104). Then, for all e e]0, 1] and for all time t e [0,r;]: 

< \ fe^'^-^^) ||£-^/:o(e,<Z^<)|||M + ||?f(V)H|M ) , 



Ci.e) 
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o Contribution of the second term. We now estimate the contribution of 

We need to clearly understand the regularity required on the pressure and the velocity. We expand 
the derivatives of the product thanks to the Leibniz formula and put aside the extremal terms in 
the summation. There exists a family of positive constant {C^.p} such that 

JtxR 



/T 

.R+M-2 



£ 

0<;3<a 



£^+*^-2 f ^2a, dev^ ^g^^ny ^^ ^y ^^3^^ 



Txl 



2 



-First, we deal with -£^+*^"^ /txm^^"' ^'^ {d'^qf^dOdy. We perform as in Inequality (129). 
For all £ e]0, 1] and for all time t e [0, T;], 



R+M-2 I 2ai ^e«e" /'Aa;?R 



£ / £■ 

Txl 



lii: 



2 



2 



(a"gf)^rf0dy 



< Jl|5«^filL==lk"^5"9f||L < ||?fil|,., . (135) 



-Then, we consider the second extremal term e^+*^ 2 j^^^^2qi goyiR Q^qR d^qf d9dy. We have 
to get a control over the family {doqf^}^. To do so, we introduce a power of e (to a cost on M). 

£^+*^-2|(e2"i a"^i«(t,.) c'egf(i,-),5"qf(i,-))| 

<£«+*^-2||9egf(t,-)|Loo |(£"^a"«i«(t,.),e"^5"§f(i,.))| , 



^fl+A/-3 / 2 2 



< 



^R+M-3 



|9?(i, oiL^.- f ii«f (i, oiiL^+s + Kit, oii]^,.. 



2 

Then we use the equivalence of norms (5) together with assumption M > 7/2, 

^fl+M-2|^^2., a«^lfl(i^.) a,5-f(i,.),9"5'f(i,.))| 

< 



^R+AI-7/2 



"■"(1,.) V -"(i.e)/ 



<(||z;f(i,-)||J,„+. + ||gf(i,-)|||Mj . (136) 

-Finally, there remains to get a bound for the sum appearing in the Equation (134). The Gagliardo- 
Niremberg's estimate only provides a control in terms of the H," ^ -norm of the pressure. Here 
one can notice that < f3 < a and so |a — /3| < to + 2. Thus d°'~^v^ is bounded in L°° . For all 



52 



£ e]0, 1] and for all time t G [0, T*], 



'e llL2 



J(l+||^f(t,.)||i„+.)||9f(i,-)||ii.| (137) 



< 

- 2 

-Finally plugging estimates (135), (136) and (137) into (134) we get for all e g]0, 1] and for all 
time te[0,T*], 

l5f(i,-)||ii=i +||«f(i,-)||i™+4 ■ 



\o<;3<Q y 



o Contribution of the fourth term. To estimate <^Ce^+*^ ^e"i9" (q^ dgvl^^ , e"i9"q^), we first 
apply the Cauchy-Schwarz inequality, 

The control over d" {q^ dgvl^^ could be done thanks to the Gagliardo-Niremberg inequality. How- 
ever it only gives a bound which depends on the -ffj? \ -norm (respectively 7Jl"l"'""'^-norm) of the 
pressure (respectively velocity). We have to be more accurate. We compute this contribution thanks 
to the Leibniz formula: 

||^gfl,+M-2^ai5a ^gR Q^ylR^ ^^ < Ce'''+"-^ ^ C„,^ \\e^' d^qf d'^-^dgvm^, . 

0</3<a 

Then we study the competition between the regularity over the pressure q^ and over the velocity v^. 
To do so, we cut the sum in two parts depending on the length of (3. When |/3| is small enough, d^q^ 
is bounded in L°° whereas when |/3| is large it is d°'~'^dgvl^ which is bounded in L°° . Therefore 
we decompose the sum into 

||^-i?.+ A/-2^Qi cia /'-R c, -IflA II ^ ^^fl+M-2 V^ ^ lU"! /I/S;?^ ,qa-/9;^ ^l-R-ll 

\\Ue £ a [q^ OgV^ j||^2 S ^£ 2^ ^a,l3\\£ O q^ O '^OgV^ \\ ^2 

O<0<afi<\0\<l 



+ Ce"+^'-^ Y. ^"^/5 Ik"' ^^^" d'^-^devl'' 






0<;9<q,2<|,9|<|q| 

•^ 
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We start by estimating C, 



Ps, 

(-^I3s < Oe 2^ ^a,p\\£ c* y£ IIloc lie "^Oev^ 11^2 

0</3<Q,0<|,3|<l 



S oe 2^ Oa^^ II g^ II 1,. 



Iw^^ll 

o<,3<Q,o<|,g|<i 



< p fl+M-2-1/2 V^ p fl ||n-R|| ll«l^ll 

0<;g<Q,0<|/3|<l '^'^' 

< 2c |e^+^^-5/2| Yl Ca,HhfL™+4- (138) 

O<0<a,O<|0|<l 

Presently, we get a bound the second term Cpi by 

1-^0/ < L^e 2^ Ca^p ||£ C^^e 11^2 ||C OgV^ ||^^ , 

0<;9<Q,2<|/3|<|a| 



< Ce«+*^-2 o„ « llC' |U,|.| llf 



<-a,/3 ||g, ll^ioi 

0</3<a,2<|/3|<|a| 






^ r^^R+M-2 V^ ^ 11 Rll ||~R|| /-,oQ\ 

< '-'^ 2_^ <-'a,/3 IPe ||//™+4 llSe ||^l°l • (139) 



0<;3<a,2<|/3|<|a| '^'^' 



Joining the two estimates (138) and (139) together we deduce that for all e G]0, 1] and for all 
time te [0,T*]: 



2 



Tl-RlP 



^a<0<a,O<\0\<l 



y \0<^<Q,2<|/3|<|a| / 



2 



This achieves the proof. D 

Finally, Lemma 3.10 is a simple corollary of Lemmas 3. 11-. . .-3.14. 

3.2.3 Proof of Proposition 3.2 

In this subsection we prove Proposition 3.2. We prove, by induction on the size J that prop- 
erty Q{J) defined page 46 is satisfied for J e JO, rn + 3] . Since the proof for the initialization of the 
induction (Q(0) is true) and the increment {Q{J) =^ Q(J-t-l)), we only the proof of the increment. 
One aspect of the proof is to give meaning to the a priori Inequality (126). 

Proof of Proposition 3.2. We assume Q( J) is true, for some J G JO, tti + 2]. 
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o Apply Lemma 3.10, there exist three positive constants Cj^i, Cj^^ and Cj_^i such that for 
all e e]0, 1] and for all time t e [0, T;]: 

at||§f(i,-)||L.+i(TxM) <'/'/^'W||9f(i,-)|||..+i +^i+\t), (140) 

with for all s e]0, 1] and for all time t e [0, T*]: 

On the one hand, since v^ is in L^HJ^I^^ (see Proposition 3.1), the family {V'/^^}ee]o,ed] is 



o 
bounded in L^ 



Vee]0,£d], Vte[0,T,*], / V/+i(s)ds<l. 

Jo 

On the other hand for all e e]0, e^], applying the assumption of induction Q{J) then Proposition 3.1, 
the function (p'^^^ is bounded in L^([0,t]) for any t € [0,t*]: 

(£2-". + ||^;f(s,.)||^„+.) ds<s'^-t. 



Jo 



We can integrate the Inequality (140) with respect to the time and apply the Gronwall Lemma. 

For all e G]0,£d] and for all time t e [0,r^*]: 



gf(i,-)^.+. < sup (^r'(s))exp / rs^'is)ds < te' 
(1-) se[Q.t] \Jo / 

It proves the induction. D 
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